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1.' INTRODUCTION 


The Ames Research Center of NASA is engaged in an effort 
to enhance its understanding of the behavior of turbulent 
boundary layers and to improve its ability to predict their 
effect on bodies in supersonic flow. As a part of this effort, 
it has entered into a contract with A.R.A.P. Under this con- 
tract, A.R.A.P. has delivered to ARC a computer program, GYC, 
which is capable of computing the properties of a compressible 
turbulent boundary layer on a rotating axisymmetric cone- 
cylinder body, according to the principles of invariant 
modeling. The program, which was delivered at the beginning 
of 1976, is to be extended during the current year to include' 
the calculation of the turbulence scale by a differential 
equation. 

In the meantime, GYC has been in operation on the ARC 
CDC-7600 computer and has undergone several corrections and 
improvements as a result of the experience gained with it 
there and at A.R.A.P. However, so far there has been lacking 
a comprehensive document giving the theoretical basis for the 
program and the method of implementation, as well as infor- 
mation on its operation. This paper is intended to fulfill 
this need, but in a tentative way. When an extended version 
of the program is delivered in a few months, it is hoped that 
a report describing it — an improved version of this — will 
be included. The anticipated improvements will be generated 
largely by the reactions to this paper. 

Section 2 is a concise review of the means used to derive 
the equations. The models used are listed but not justified. 

^ Some background on the models can be found in Reference 1 for 
the models that apply to Incompressible flows, and in Refer- 
ence 2 for others . 
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In Section 3 , the numerical method used to solve the 
equations is set forth. Some extra computations are described 
in Section 4, while the means for controlling the finite- 
difference mesh are reported in Section 5. 

Sections 2 through 5 are essentially free of program 
considerations, which are covered in Sections 6 and 7* Section 
6 describes GYC as a whole and then considers in more or less 
detail the main sections of the program. The part concerned 
with the numerical integration of the equations is given the 
most attention. External aspects of the program operation 
are described in Section 7- The paper is brought to a close 
with a few remarks in Section 8. 

The full equations under consideration are set forth in 
appendices. An extensive glossary of Fortran names and Fortran- 
like names is included as Appendix F. 
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2 . ANALYSIS 

The basic equations used in this study are the following: 
Continuity 


Momentum 


where 


pu. 


+ pu u 


i,Jl 



+ 


T 


£ 

i,£ 


T'! = y*5^'u^ 

X \1 ^1/ X3 


Enthalpy Version of Energy 

ph^ + pu ' Pt + u p.d + 

The thermodynamic relations of a "calorically perfect" gas are 
used 

p = pRT 

h = c T 
P 

where R and Cp are constants . 

Writing the dependent variables in these equations as the 
sum of mean and fluctuating parts ( p = p + p' , for example), 
we can deduce equations for the mean quantities and for the 
second-order correlations of the fluctuations by protracted 
manipulations. Two assumptions are made immediately: 

1) Fourth-order correlations are neglected. 

2) Third-order correlations involving y , y* 3 or k 
are neglected. 
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The resulting equations for a steady mean flow are displayed 
in Appendix A. 

Fluctuations in the transport parameters are handled by 
assuming 


y' = UrpT* 


= U|T' 

u*;i - 

k’ = k^T’ 

k' = kjT'. 


where the subscript T denotes a derivative with respect to 
temperature evaluated at T . 

The models used to close the equations are given in 
Appendix B. As in previous A.R.A.P, models, these are 
invariant under coordinate and Galilean transformations and, 
of course, are dimensionally consistent. The goal of sim- 
plicity has been sacrificed to some extent in favor of 
generality. An obvious feature of these models is that many 
of them are expressed in terms of quantities that are them- 
selves modeled. This results in some of the parameters 
(coefficients) appearing only in products with other parameters 
in the applications. (This is exemplified by the models for 
p ' h ' • and p ’ p ' • where the substitution has already been 

5 -L 

carried out.) Thus some of the parameters are redundant, but 
they are left that way for the sake of flexibility in case it 
is found necessary to modify them in the future. 

Three- or four-character names, suitable for direct use 
in Fortran, have been chosen for the modeling parameters 
except for a , b , and B , in the velocity-dissipation 
model, which are left in that form out of nostalgia. 

It will be seen later that all of the terms containing 
models for correlations involving the divergence of the velocity 
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fluctuations, u'^ , drop out. For this reason, the following 
parameters do not appear in the final equations: BDV, VBS, ¥MD, 

WWU2, WWDl, WWD2, WWH2, and ¥WR2 . In addition, terms involving 
PMH, PMR, PTHM, and PTUM do not occur in the final set. Of the 
remaining 38 independent parameters, about a quarter are set 
to zero until better information on them is available, and the 
majority of the others are evaluated by analogy. Thus 0.1 has 
been found a good value for VUU in constant density tests; VRU, 
VUH, VRR, VRH, and VHH don’t occur in constant -density modeling 
but, by analogy, 0.1 is used for all of them, until better infer 
mation is available . 

So far the equations have been written using mechanical 
units for heat quantities. In order to introduce additional 
flexibility, a conversion factor is introduced, allowing arbi- 
trary units for heat. This factor is 

GMOMS = 1/J 

where J is the mechanical equivalent of heat. For example, 
if English units are used with heat measured in Btu, GMOMS = 
0,001285 Btu/ft-lb (= 1/(778 ft-lb/Btu)). If mechanical units 
are used for heat, GMOMS = 1 . It is often desirable to make 
the equations nondimensional . If this is done using u^, and 
hp as reference values for velocity and enthalpy per unit mass, 
the nondimensional version of GMOMS is 

2 2 

u u « 

GMOMS,, ^ = GMOMS ~ = = (y - 1)M 

ND h Jh ^ ' r 

r r 

(Hence the name GMOMS - Gamma Minus One M^ Squared.) 

The equations with the models substituted are exhibited in 
Appendix C. Here the notation is that devised for input to 
TENSR, A.R.A.P.'s software for the expansion of tensor equations 
This notation is like Fortran in several respects. Thus, 
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variable names may consist of more than one character; multipli- 
cation is always made explicit with a * ; multiplication and 
division proceed from left to right so that A/B/C*D is to be 
evaluated as ((A/B)/C)kD = (A«D)/(B*C) . On the other hand, 
superscripts and subscripts are represented in a manner com- 
pletely different from the Fortran repre'sentation of subscripts. 
Superscripts are preceded by " ; subscripts by ’ . Super- 
scripts and subscripts consist of a single letter or digit only; 
two or more together represent separate indices. The metric 
tensor is denoted by @ ; the Kronecker delta by & . Covar- 
iant differentiation is denoted by a , preceding an index; 
and "contravariant differentiation" (that is, covariant differ- 
entiation with the index raised) by I preceding an index. 

All these forms are illustrated in Appendix E. 

The Glossary (Appendix P) gives the meaning of the variable 
names used. Thus, for RS, US, and -RU in the first equation of 
Appendix C, we find p , u , and p 'u* respectively. With 
the help of Appendix E, we see that the first equation of Appen- 
dix C is a transcription of the first equation in Appendix A, 
as it should be since there are no terms needing modeling. 
Proceeding thus step by step, the formidable expressions of 
Appendix C can gradually be seen to make sense. 

A cross-section of the axisymmetric coordinate system is 
shown in the sketch on the following page. 
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= y = pQ ^ = ({) = the azimuthal angle. The half angle 

of the cone, 0^ , and the "origin radius," Tq , are 
parameters of the system. The coordinates x and y are 
related to the cylindrical coordinates ^ = SR and r = PR 

by 

5 = cos 0 qX - sin QqY 
r = r^ + sin 6j,x + cos Q^y 

It is understood that y = 0 always represents the cone sur-- 
face. The location of the origin in the x direction is quite 
arbitrary and hence, in general, the value of rQ is arbitrary. 
The exception is the limiting case of a cylindrical surface 
(0„ =0 or e„ = 180°) for which r^ represents the radius 
of the cylinder. 

Note that for 8^ in the second and third quadrants, the 
flowfield (for which y > 0 ) is on the inside of a conical 
surface . 
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It is convenient to have the equations apply also for a 
Cartesian coordinate system (x,yjz). To accomplish this, 
another parameter is defined by 


Now let 


DIMPL 


fl. axisymmetric 
t^O. Cartesian 


S = DIMPL* sin 6^^ 

C = DIMPL^t cos 0^ 

R = DIMPL*rQ + Sx + Cy + 1. - DIMPL 


Then the metric tensor can be written 


10 0 


gij = 1 0 1 0 

‘ 0 0 R" 


The nonvanishing Christoffel symbols of the second kind are 

$33 = -SR 

$33 = -CR 

$31 = $13 = S/R 

$ 3 (t 3 — r* /"D 

32 ^23 o/n 


These apply for both Cartesian coordinates and the axisymmetric 
coordinates of the sketch. 


8 



Note that 

{ r axl symmetric 

1. Cartesian 

so that it not only has different meaning but even different 
dimensions in the two cases. 

The next stage of the analysis was carried out by TENSR, 
A.R.A.P.'s software for the expansion of tensor expressions. 

The following information was supplied to TENSR: 

1) The values of the metric tensor and the Christoff el 
symbols of the second kind in terms of R , S , and 
C , as given above . 

2) The information that terms involving derivatives with 
respect to x^ (^) are zero^ since the flow is 
assumed to be axisymmetric . 

3) The substitutions required so the output involving 
the mean velocity could be expressed in terms of the 
physical components, u , v , and w (or US , 

VS , and WS ) , Instead of the covariant and contra- 
variant components, and similarly for correlations 
involving u ’ . 

4) The equations of Appendix C. 

The output for the whole set of l8 equations consisted of 
98^3 terms. This number is somewhat deceptive, since TENSR is 
naive algebraically and does not combine like terms. 

These equations were submitted to more A.R.A.P. software 
known as TATTR. The first stage of TATTR does combine terms 
and reduced the l8 equations to a total of only (!) 43^8 terms. 

The main function of TATTR is to drop terms according to 
order-of -magnitude assumptions. Its operation in this case can 
be described as follows. 
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The equations were interpreted to be written in non- 
dimensional form. The (nondimens ional) boundary-layer 
thickness was assumed to have a small value, 6 . The mean 
quantities and the radius, R , were assumed to be of order 
one except for the normal component of the velocity, v (VS) , 
which was taken to be of order 6 . Derivatives with respect 
to X and y were assumed to be of order one and 6 , 

respectively. The transport coefficients p , u* , and k 
were assumed to be of order 6^ . All of the above corresponds 
to standard boundary layer' assumptions. In addition, correla- 
tions involving a velocity fluctuation were taken to be of 

1/2 1/1 
order 6 ' and those Involving p’ or h' of order 

for each such appearance. (For example, u’u’ , u’v* , p *w’ , 

and h'h’ ( UU , UV , RW , and HH ) were assigned orders 

6 , 6 , and , respectively.) The scale length 

A was assumed to be of order 5 ^^^ . Finally, GMOMS , which, 

p 

as has been seen, is (y - 1) when the equations are inter- 
preted nondimenslonally , was taken to be of order , 

All quantities not mentioned. Including the modeling- parameters , 
were assumed to be of order one. These ratios v^ere arrived at 
after some trial and error, using Cartesian coordinates. They 
were chosen so as to insure the retention of terms known to be 
Important, while allowing the elimination of terms known to be 
unimportant . 

TATTR computed the exponent, T , representing the order 

rp 

of magnitude 6 , for each term of each equation, given the 

information in the paragraph above concerning the factors that 
can appear. For example, the continuity equation was analyzed 
as follows (where :1 and :2 stand for partial differentiation 
with respect to x and to x ) . 
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T 

(RS*US):1 0 

+ (RS«VS):2 0 

+ S»RS/R*US 0 

+ C*RS/R*VS 1 

+ (RU):1 5/6 

+ (RV):2 -1/6 

+ S/R«RU 5/6 

+ C/R«RV 5/6 

= 0 


The final stage of -TATTR dropped terms for which the value 
of T is greater than a maximum appropriate for the particular 
equation. In the case of the continuity equation, the cutoff 
value is zero, so the final form is 

(RS*US):1 + (RS*VS):2 + S«RS/R»US + (RV):2 = 0 
or 

(pU)^ + (PV)y + fpU + (P'V’ )y = 0 

The whole set of equations was reduced from 43^8 terms to 
529 terms by the final stage of TATTR. The second component 
(normal to the surface) of the momentum equation reduced to 

T 

-C»RS/R«WS«WS 0 

+ (RS*W):2 0 

= -(PS) :2 -1 

The values of T indicate that, according to TATTR, the 

remaining terms are seriously out of balance. But it is well 
knovm in boundary-layer theory, though not to TATTR, that the 
derivative of pressure in the normal direction is not large 
like other normal derivatives; (PS): 2 is really of order one 

( 5 *^) , not 6"^ . Therefore, the terms containing (PS):2 in 
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the other equations were dropped manually if the T value for 
that term was within one of the cutoff value for the equation. 

In this way, 63 terms were dropped leaving 466 in the 18 equa- 
tions. Two of the final equations have been presented above. 

The others, which have been slightly edited manually, are pre- 
sented in Appendix D. (The reason for the organization of the 
equations into "passes” is given in the next section.) 

As has been indicated, the criteria used to' drop terms 
were arrived at pragmatically. This can be defended on the 
basis that (1) the process is more systematic than that hereto- 
fore used for modeled equations for compressible turbulent flow 
and (2) the results are not sensitive to the details of the 
assumptions used — e-g., it makes little difference whether p' 
is assumed to be of order 6^'^^ or order 6^^^^ . Nevertheless, 
it would be valuable to find a rational basis for completely 
specifying the criteria. Such a basis should take into account 
the fact that the balance of terms differs markedly in different 
parts of the' boundary layer. 

Some features of the final set of equations should be 
noted. (1) All terms involving models for u^£ correlations, 
as noted above, have dropped out. This is a result, of course, 
not only of the criteria used, but also >of the particular models 
used for such correlations. (2) All terms involving the second 
coefficient of viscosity, y* , have dropped out. (For the 
most part, these are terms that also Involve .) (3) If S 

and C are both zero (as they are if DIMPL =0), the equations 
consisting of the remaining terms are each homogeneous in R . 
This fact (which is true of the equations before terms are 
dropped) shows that the change in dimensionality of R as 
DIMFL changes from one to zero does not violate any principles. 

A new equation 


BLAM = GYCSZ 
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has been added to the set and included in Appendix D. Its 
presence here is merely symbolic. BLAM (A) is evaluated as 
a known function of y and the gross features of the flow, 
quantified by input parameters (the scale parameters). It is 
included here and in the program as a pseudo-unknown, satis- 
fying a "differential" equation of zero order which becomes 
a "difference" equation of zero order. The reason for this 
treatment is that it is intended to introduce a true differ- 
ential equation for A and incorporate its solution in the 
program in the near future. The program modification will be 
much simplified by having A already treated as an unknown. 

The sixteen differential equations of Appendix D are 
parabolic in the sense that they have first-order derivatives 
with respect to x and second-order derivatives (as well as 
first-order) with respect to y . Thus x is the time-like 
variable and y is the space-like variable. The dependent 
variables determined by the equations are referred to as para- 
bolic variables. 

The boundary conditions needed for this set consist of 
the initial conditions, the wall conditions, and the free- 
stream conditions. The Initial conditions consist of profiles 
of the parabolic variables, that is, their values as functions 
of y at the x station where computation is to start. The 
wall conditions on all the turbulence correlations and on u 
are that they are zero there. The wall condition on w is 
that it is equal to the lateral velocity of the cone surface, 
that is, the product of the angular velocity and the local 
cone radius- The wall condition on h is either its value at 
the wall or the value of its gradient (a measure of the heat 
flux) at the wall. Either of these may be functions of x . 

The free-stream conditions, or "edge conditions," are that 
the parabolic variables approach given values asymptotically. 
For the turbulence correlations and for w , the given values 
are constant in x , but for u and h they may be functions 
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of X . In practice the concept "approach asymptotically" is 
handled as follows. At the largest value of y currently used 
in the computation, the boundary condition ^ = 0 , where f 
stands for any of the unknowns, is used in solving the equations. 
The resulting value of f at that y is compared with the 
given value and if they agree within a specified tolerance, the 
given value is substituted for the calculated value (for the 
sake of uniformity and neatness) and the solution continues. 

If the tolerance is not satisfied, for any of the parabolic 
variables, the maximum value of y is increased by adding 
another point, and the process is repeated. 

These boundary conditions are not the most general that 
can be posed for the system of equations; rather they are the 
conditions provided for in the GYC program. 
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3. NUMERICAL METHOD 


A finite-difference mesh is introduced in the x , y 
plane as shown in the sketch. 



The time-like variable, x , is indexed by n and the space- 
like variable, y , by j . The spacing is not uniform in 
either direction; furthermore, the spacing in the y direction 
is adjusted as the solution proceeds in the x direction, as 
described in Section 5. 

At the beginning of each x step, the values of the 
dependent variables are known at x = for all y^ . 

Derivatives with respect to x are approximated by 


where 




15 



some additional notation 




For. derivatives with respect to y , 
is convenient: 


h+ = yj+i - yj 


- yj-i 


' yj+i - yj-1. 


H 




h, 


H ' 


h_h^. 


0^.(8) 


n-1 , n-1 

^ 3+1 


G (g) 


n-1 . n-1 

_ S i S ,_1 

h_h^ 


Then we have 


"y = - ^j) " «-("” - ‘'d-i) 

K)y “ °+<8> (^J+l - ‘■j) - - ^3-l) 

If g is constant, equation (3) can be written 

(^^1 - ^f) - h4r(1 - 


■yy 


( 2 ) 

(3) 

(4) 
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One way of arriving at equations (2) and ( 4 ) is by passing a 
parabola through the three points ^ ) ’ 

(^j+1’ ^j+l) evaluating its derivatives at yj . Equation 
( 3 ) is an obvious generalization of equation ( 4 ) except that g 
is evaluated at instead of at . This is done to 

linearize the numerical procedure since g , in general, depends 
on the unknowns. The other nonlinear terms in Appendix D are 
handled similarly; that is, in any product, at most one factor 
is evaluated at x„ while the rest of the term is evaluated at 


'n 


X 


n-1 


This artifice has been successfully used in similar 


equations for several years at A.R.A.P. Of course, it is neces- 
sary to keep Ax small enough- so that the changes in the depend- 
ent variables are small in each step. The linearized system is 
in the form referred to in the literature as the (complete, or 
fully) implicit method for the numerical solution of parabolic 
differential equations. 


For terms of the form pvfy , equation (2) is optionally 
replaced by an, upwind differencing formulation, described in 
Note 4 of the Glossary (Appendix P) . 


The equations are coupled, but the finite-difference 
versions can be decoupled to any extent desired by the same 
device that is used to make them linear. That is, terms or 
factors which are unknowns if evaluated at x^ become known if 
evaluated at . Further, if some unknowns are solved for 

before others, they become known at both and x^^ for 

the later solutions.. If m equations are solved at once, at 
each point of the flowfield a matrix of size m • m must be 
inverted and two matrix multiplications of m • m matrices 
must be performed. Since these procedures require of the order 
of m^ operations each, there is strong motivation for keeping 
m down. On the other hand, experience at A.R.A.P. has shown 
that some equations are better -solved together in order to have 
stability and well behaved solutions with reasonable step sizes. 
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For GYC, these considerations have led to the division of 
the solution process, for each step in x , into six passes. 

In the first pass, w and v’w' are determined so m = 2 . 

In the rest of the passes m = 3 : u , u* v’ , and u*w* are 
determined in the second; u 'u ' , v’v' , and w'w ' in the 
third; h , h'v' , and h'h' in the fourth; p'u’ , p’v' , 
and p'w’ in the fifth; and h’u' , h ' w' , and A in the 
last. 

In each pass, the m unknowns are thought of as forming 
a vector, (j)^ , of m elements at each yj (for x = Xj^ ) . 
Thus in the first pass 



in the second 



etc . 

* 

Then the finite-difference equations can be written 






(5) 


where Aj , Bj , and Cj are square matrices of order m 

and d,- is a vector of dimension m . The elements are known 
J 

functions of the dependent variables at (and in a few 

instances at x^^ for variables determined on previous passes), 
the independent variables, and, through equations (l)-(M), the 
mesh spacing. 
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Equation (5) applies for 2 < j < J - 1 where J = JTOP 
is the index for the largest y in use and J = 1 corresponds 
to the wall (y-|_ =0) . The finite-difference version of the 
edge boundary condition, that the derivative with respect to y 
be zero for each variable, can be expressed as 

'^J-1 ~ 


The wall boundary conditions can be written 


<i>i 


d 


1 


(7) 


except when h is involved and its gradient at the wall is 
specified. 

When the gradient of h at the wall is specified, the 
boundary condition is handled this way. It is required that 
there be no blowing so both components of the mean velocity, 
as well as all correlations, are zero at the wall. Then, on 
the wall the differential equation for h (the first equation 
of the fourth pass) reduces to 


0 


U 


-2 


+ (Rw) 

y 



+ 



y 


assuming mechanical units so GMOMS = 1 , The condition on 
the gradient is written 



= g 


where g is the specified value. These two equations are 
written in finite difference form using three points, the wall, 
y^ = 0 , the first point off the wall, y2 a and a "phantom" 
point, yg , inside the wall. (Extrapolation is used to 
evaluate "known" quantities at the phantom point.) Thus there 
are two equations involving the unknowns hg , h^ , h^ . 
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Eliminating between these equations leaves one equation 

which can be written 


ahj + 3hg = T (8) 

where a , 3 , and y depend on the known variables and the 

mesh spacing at the wall. 

In any case, we conclude from equati'ons (7) and (8) that 
the wall boundary conditions for a pass can be written 

® 1^1 ^ 1*^2 " ^1 
< 

where is diagonal and Cj is zero except in the fourth 

pass when the gradient of h is specified, in which case the 
upper left element of is the 3 of equation (8) and the 

rest of C]_ is zero. 

Similarly, equation (6) can be written 

Vj-1 “ '^J 


where Aj is the unit matrix, Bj is the negative of the unit 
matrix, and dj is zero. 

Equations (5)j (9)j and (10) can be written 
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or 


M$ = D 

where M is a tridiagonal matrix of order J , each element 
of which is an m • m matrix, and $ and D are vectors of 
dimension J , each element of which is a vector of dimension 

m . 

This set is solved by a standard algorithm for tridiagonal 
systems. The process consists of two parts called the upsweep 
and the downsweep. In the upsweep the values of auxiliary 
matrices, Fj , and auxiliary vectors, , are computed by 

d^ 

"■j = J = 2, 3. .... J-1 (11) 

° ^('*3 ■ *J*d-i) i ' 

In the downsweep the solution is obtained by 

“ c()j — ” J“1 , J—2 , 1 (13) 

( A-j_ and Cj do not appear in the system. If we define 

= Cj = 0 , equations (11) and (12) can be used for J =1, 2, 
..., J ; equation (13) can be used for j = J, J-1, -..j 1 ", 
and the special equations for 0=1 and j = J are not ' 
needed. ) 
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As mentioned in the previous section, the edge boundary 
condition is really asymptotic. The imposition of equation (6) 
or (10) at yj is only tentative. At the end of the upsweep, 
the values of tjjj = are compared with the known free-stream 
values. If they are not equal within specified tolerances, , J 
is increased by one; the already known solutions are suitably 
extrapolated; Aj , Bj , , and d^ for j = Jq , the 

old value of J , are recomputed using expressions valid for 
j < J , the upsweep is redone for j = Jq and done for j = J ; 
and the test is repeated. When the tolerances are satisfied, 
the given free-stream values are substituted for the computed 
values at J (from which they differ by less than the toler- 
ances) and the downsweep is performed. 

Consider a regime where w = 0 , the turbulence is negli- 
gible, and y derivatives are negligible. The equations for 
the means reduce to 

puUx = -Px 

puh^ = up^ (15) 

(taking GMOMS = 1 ) . If the edge conditions on u , h , and 
p satisfy these relations ( p being determined by the equation 
of state), this handling of the free-stream boundary -conditions 
works well. Otherwise, an unreasonable number of points may be 
added, satisfying finite-difference versions of equations (l4) 
and (15) at the added points while trying to match conflicting 
given edge values. In this case, it is probably not right to 
assume y derivatives are negligible; but if they are not, it 
isn’t clear how the edge of the calculation field should be 
determined . 

If there is free-stream turbulence, another situation 
arises. Consider the simplest case with w = 0 , the other 
mean quantities constant, y derivatives negligible, and 
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u‘^u*. = 6^q^/3 

J J 

2 

Then it is found that q is governed by 



That is, dissipation will make free-stream turbulence decrease 
downstream unless A -»■ <» , according to the equations being 
used. There ia no provision in the program as it stands for 
the turbulence-level edge condition to vary with x so exces- 
sive points will be added if free-stream turbulence is specified 
unless A is made to increase for large y . The scale param- 
eters which control the calculation of A allow for such 
increase. 

Excessive adding of points at the edge can always be 
avoided by increasing the tolerance used to check the matching 
of edge conditions for the variables giving trouble. This may 
give results of questionable validity. 

There are three dependent variables that are not deter- 
mined from parabolic equations . These are p , p , and v . 
The mean pressure is calculated at x = x^^ during the down- 
sweep of the first pass from the normal momentum equation 
displayed in Section 2. The integration starts with a given 
edge value of p (which may be a function of x ) and proceeds 
inward to the wall. The mean density is calculated at x = Xj^ 
during the downsweep of the fourth pass from the mean equation 
of state, which can be written 


R(h - h’hVh) 

The normal component of the mean velocity is calculated at 
X = Xj^ after the completion of all the passes, using the 
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continuity equation which is also displayed in Section 2. The 
integration starts with a given wall value of v (which may 
be a function of x ) and proceeds outward to the edge. Trape- 
zoidal integration is used for both p and v . 

Another quantity, namely u„ , is grouped with p , p , 
and V . Sinc'e its value for the current step is not available 
in the first pass, its value saved from the previous step is 
used there. The contribution of the terms in which it appears 
is expected to be small, so the error involved in lagging it 
this way should also be small. The value from the current step 
is set aside for use in the next step as part of the downsweep 
of the second pass. 
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H. AUXILIARY COMPUTATIONS 


Some functionals of the dependent variables are calcu- 
lated at the completion of each x step (as well as during the 
initialization stage). These calculations serve three main 
purposes: (1) they provide output which helps the user interpret 
the results of the program, (2) they provide data used in deter- 
mining the finite-difference mesh for the next step, and (3) 
they provide data on which the calculation of the local value 
of A is based. In the last case, the process is not really 
auxiliary since A is basic to the modeling and strongly influ- 
ences the course of the solution. 

Many of the calculations are repeated for each of the 
parabolic variables. To avoid long descriptions, such proce- 
dures are set forth here as if there were only one Independent 
variable, named F . 

These quantities include: (1) the "max" which is the maxi- 
mum over yj of ]Fj at ^ value of y at which 

the max occurs, (3) the "global max" which is the maximum over 
yj and x^ (up to the current x^ ) of |P| , (4) the "max 

change" which is the maximum over y^ of 

(5) the integral (according to the trapezoidal rule) 

nj 

1(F) = I iFjdy 

J o 

(6) the "integral spread" which is 1(F) divided by the max of 
F , and (7) the "percent spread" which is the value of y for 
which F is a designated fraction of its max. 
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In addition, , which is the value of y for -which 

u is 99 percent of its free-stream value, is found. A typ- 
ical length, , is defined as a linear 'combination of 599 

and the percent spread of one of the parabolic variables. The 
coefficients are inputs as is the choice of which percent 
spread to use. 

The calculation of A , the turbulence scale, is given by 



^i = 

min ^ 




A = \ 

h 


y - 'S99 



+ 

- S9) 

y > <599 

The 

quantities c , 


d^ , and Aq 

are input parameters. 

For 

normal values of 

these 

constants, A 

equals Aq at the 


wall and increases with slope dj^ until it reaches the value 
cL^ . It stays at that value the rest of the way out if d^ 
is zero. Otherwise, it starts increasing again at y = 6 ^^ 
with slope d^ . The parameter Aq is made greater than zero 
to simulate wall roughness. The parameter dQ is made- greater 
than zero to allow for free-stream turbulence without dissi- 
pation there, as mentioned in Section 3. Otherwise,, both A^ 
and d^ are normally zero . 

The rest of the auxiliary calculations are done to provide- 
informative output only. They have no effect on the course of 
the program. The output labels and definitions are: 


momentum thickness 


THETA 




dy 
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displacement thickness 


DELTA* 




shape factor 


SHAPE = DELTA*/THETA 


kinematic displacement thickness 


wall shear 




TAUW 



y=o 


local skin friction coefficient 


CF = TAUW/(|pgUg) 


wall heat flux 


QWALL 


Cp sy 




2k 


total Stanton number 


STANT = QWALL/ 


^e e 




Reynolds numbers 


RETHETA = THETA 


P ^ 
^e e 


REXW = XW 




laminar stress 


TAULAM = 


- 3u 


total stress 


TAUTOT = 


- 3u 


pu'v' 


total heat flux 

QTOT = 

The integrals are done by the trapezoidal rule. The derivatives 
are evaluated by equation (2) of Section 3 except at the wall 
where a two-point formulation is used. In this regard it should 
be pointed out that the u term in the formula for QWALL is 
superfluous as a derivative (since u = 0 at the wall, the 
derivative of its square is zero there), but in finite-difference 
terms the form given can be shown to be more accurate, at least 
if w = 0 and v = 0 at the wall. 

None of the calculations in this group take into account 
the possible effects of nonzero w . Modifications or addi- 
tional calculations to show such effects could easily be intro- 
duced into the program. 
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5. CONTROL OF THE MESH 


In order to provide reasonable accuracy without excessive 
computing; the mesh is modified by the program as the run pro- 
gresses. Parameters supplied by the user allow him to control 
the balance between accuracy and computing time. 

The value of Ax for the first step of a run is supplied 
as an input. ’Thereafter, a new Ax is computed before each x 
step. To do this, tentative values of the step size are com- 
puted for each of the 17 parabolic variables (that is, the 
variables determined by parabolic equations), and the smallest 
of these is selected. Each tentative value is determined by 
multiplying the previous Ax by the ratio of a criterion estab- 
lished for the particular variable to the max change of that 
variable for the previous step. The criterion is a nominal 
change that can be specified either as an absolute quantity or 
as some factor times the global max of the variable with an 
optional lower limit . The new Ax , the minimum of the tenta- 
tive values, is then modified if necessary so that it doesn't 
exceed either a specified factor times the old Ax or a speci- 
fied . On the other hand, if the new Ax is less than 

a specified fraction of the old, "backup" occurs.’ That is, some 
change generated by the previous step is considered unacceptably 
large so the results of the solution for conditions at Xj^ 
just computed are discarded, and a smaller step from the con- 
ditions at is tried, using the new smaller Ax . A 

minimum Ax may also be specified. If Ax becomes less than 
the minimum, the run is halted. 

A further complication results from the existence of so- 
called time breaks. These are values of x at which the user 
requires output. In order to avoid an abnormally small Ax as 
the last step to such a point, the program looks several steps 
ahead and modifies the Ax found by the process described above, 
if need be. 
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The initial distribution of points in the y (normal) 
direction is part of the input of the initial conditions. 

Before the first and before each succeeding x step^ a complex 
process is used to monitor and, if indicated, to adjust the 
spacing of points in the y direction. The basic idea is to 
examine each set of three adjacent points, yj-i , yj , yj+1 • 
If, for all the dependent variables, the values for those 
points fall almost on a straight line, the mid-point can be 
dropped. On the other hand, if, for at least one of the depend- 
ent variables, there is an excessive departure from a straight 
line, a point should be inserted in one of the subintervals 
yj_l j yj 02? yj , yj+]_ , or in each. (These notions are 
made quantitative by comparing a nondimenslonal measure of the 
departure from a straight line, for each dependent variable, 
with input parameters. The measure is the absolute value of 
the second derivative, as given by equation ( 4 ), divided by the 
global max and multiplied by the square of the length of the 
appropriate interval.) However, it is also required that the 
ratio of the lengths of adjacent intervals be less than a speci- 
fied value. This requirement can prevent a point from being 
dropped, that otherwise would be, but does not prevent an inser- 
tion. Instead, insertions are made in other intervals until the 
ratio requirement is satisfied everywhere. 

There are also means of specifying maximum and minimum 
interval lengths. The maximum prevents the dropping- of a point 
which would cause it to be exceeded. If the insertion of a 
point creates an interval length smaller than the minimum, a 
warning message is printed. 

Furthermore, if the insertion of points would make J = JTOP 
bigger than its upper limit, JMAX (an input) , the user has 
the option of allowing the relaxation of the criteria governing 
the dropping and insertion of points. If the relaxation is not 
allowed or is unsuccessful in keeping j i JMAX » "the run halts. 
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When a point is inserted. It is not necessarily put at 
the mid-point of the interval. Rather it is weighted toward 
the point with the largest (in absolute value) normalized second 
derivative considering all the dependent variables. The restric- 
tion on ratio between the length^s of adjacent intervals is 
observed, of course. 

The values of the dependent variables that are assumed for 
an inserted point are obtained by linear interpolation. This 
seems like an anticlimax, but it is the only simple scheme that 
has been found here that does not lead to undesirable results 
in some cases. 
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6 . PROGRAM STRUCTURE 


The principles of invariant modeling have been applied at 
A.R.A.P. to several different types of turbulent flows that can 
be described by parabolic equations. The technique developed 
for solving those problems has been used for the GYC program] 
and, to a large extent, the actual subprograms have been carried 
over, with varying amounts of modification. 

What has been referred to as the GYC program is actually 
a main program, GYCON, and a number of subprograms. GYCON 
itself is quite simple. Basically, it consists of a call to 
GYCIN, which handles input and initialization, followed by a 
loop consisting of calls to; (1) GYCAL, which hahdles the 
auxiliary computations, the determination of the new Ax , and 
output; (2) GYCAA, which handles the distribution of the points 
yj ; (3) GYCSS, which does the actual solution for one x 
step; and (4) GYCRV, which determines v . The call to GYCRV 
is followed by a return to the call to GYCAA. This basic flow 
of the program is modified by certain special conditions; for 
example, if backup is required (see Section 5), the call to 
GYCAA is skipped. When the run is halted, GYCIN is called' 
again. Depending on the input cards supplied, this can result 
in a restart of the same run, the start of a new run, or job 
termination. 

Before describing the operation of the principal sub- 
programs mentioned above, an aspect of data handling should be 
mentioned. Some of the other programs from which this one 
evolved needed more data 'Storage during the solution process 
than was available in core (especially when the A.R.A.P. com- 
puter had only 8,192 words of 16 bits each in core). Such 
programs used external (disk) storage; and the GYC program, 
although contained entirely in core, mimics that arrangement by 
keeping a good part of the working data in a two-dimensional 
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array called STOR. The elements of STOR are never addressed 
directly; instead, parts of STOR are copied to and from other 
arrays in core, much as if STOR were an external device. 

STOR consists of 52 elements for each j , that is, for 
each normal mesh point. The first element is the value of yj 
itself. The next 17 elements are the old values of the para- 
bolic variables, that is, the values evaluated at (x^-i > yj ) 
in the order they are listed in Table 5* Then come the old 
values of p , p‘ , v , and u^ . The following 17 elements 
are used for the new values of the parabolic variables and also 
for the values of <f>j , in the notation of equation (12), 
while generating the new values. The next four elements are 
the new values of p , p , v , and u^ . The rest of the 
52 elements are used to hold the matrix fj . 

The copying or moving in and out of STOR is done by a 
subprogram, SFVMV. With parameters A, B, and N, SFVMV is 
equivalent to the simple DO loop: 

DO 8 I = 1,N 
8 B(I) = A(I) 

It exists as a subroutine because it is used so often; at A.R.A.P. 
it is implemented in firmware to speed its operation. A com- 
panion subprogram, SFVFL, fills N elements of the vector B 
with the value of the variable A , when the parameters are 
A, B, and N. 

Of the principal subprograms, GYCIN is the most straight- 
forward. It initializes certain quantities, it reads the first 
data card, it calls GYCRI to read the rest of the data cards if 
the first card Indicates that there are more, it calls GYCIC to 
do more initialization if a new run is being started, and 
finally it calls GYCPI to print the title pages of the output. 
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GYCAL first calls' GYCQD and GYCAQ which do the auxiliary 
computations described in Section 4 and which also compute the 
new Ax as described in the first part of Section 5, except 
for the modification due to time breaks which is done by a call 
to GYCTB, The bulk of GYCAL is concerned with calls to GYCOI 
and GYCOT’ if the various criteria so indicate. GYCOI produces 
a minor printout, whereas GYCOT produces a major printout 
including tables of the dependent variables as functions of y . 

GYCAA is the most complex of the subprograms. It is also 
the least changed from versions in use for other problems. Its 
function is to monitor and, if Indicated, to adjust the spacing 
of points in the y direction, as described in Section 5- 
Suffice it to say that it works as advertised. 

GYCSS, together with its subprograms, forms the heart of 
the system. It performs the upsweep and the downsweep, repeating 
the process for each of the passes described in Section 3* Since 
three-point differencing is used in the y direction, a section 
of COMMON is used for values of variables at j-1 , that is, 
variable names ending with the letters M and I (compare the 
sketches at the beginning of Section 3 and the beginning of the 
Glossary, Appendix P) ; another section of COMMON is used for 
values of variables at j , with names ending With the letters 
Z and J ; and a third section is used for values of variables 
at j+1 , with names ending with the letters P and K . 

As j is incremented at the beginning of a stage of the 
upsweep of a typical pass, the values in the Z-J section are 

shifted to the M-I section, and those of the P-K section to the 

Z-J section. Then GYCPP moves to the P-K section from STOR the 
values that are available there and calculates the values of 
others. Obviouslj, special provisions are made for j = 1 . 

The values of the elements of Aj , Bj , Cj , and dj are 

calculated by the subroutines GYCBI (for j = 1 ), GYCMC (for 

2 < j 1 J - 1 ), and GYCBO (for j = J ). After making, some 
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calculations common to all passes, GYCMC in turn calls separate 
subroutines for each pass. (Due to some confusion in the early 
stages of the development, the subprogram for the first pass is 
named GYCP2 and that for the second pass is named GYCPl.) It 
is in these pass subprograms that the actual equations being 
solved (Appendix D) make themselves felt. They were automati- 
cally generated, for the most part, by DIFFR, another A.R.A.P. 
software package . 

DIFFR took the equations of Appendix D, one pass at a time, 
along with some information about which symbols to treat as 
unknowns, ‘as known variables, or as constants, and generated 
Fortran statements evaluating the elements of Aj , , Cj , 

and dj . These statements, after some minor manual editing, 
constitute the bulk of the pass subprograms. 

Returning to the description of a stage of the upsweep, 
fj and (equations (11) and (12) of Section 3) are evalu- 

ated using the special matrix subroutines GYCMP and GYCMI. 
Finally fj and are moved to STOR, and j is incremented 

again . 

The handling of the free-stream boundary conditions has 
already been described. When they are satisfied, the downsweep 
according to equation (13) of Section 3 is trivial. It is 
slightly complicated in the first, second, and fourth passes by 
the incorporation of the calculation of p , j and p , 

respectively . 

There are several other subprograms which have not been 
mentioned here. All are Included in the Glossary (Appendix F) 
with brief descriptions. At the beginning of the listing of 
each subprogram are comments giving another view of its role. 

Three ficticious subprograms are included in the Glossary, 
and places for calls to them are suggested by comments in the 
program. These provide for the possibility that it may be 
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desirable to create files in external storage for later plotting. 
Two types of files are anticipated. One, the "running file," 
would be of quantities that depend on x only, such as skin 
friction and momentum thickness. The other, the "profile file," 
would be of quantities that are functions of both x and y , 
such as any of the dependent variables. The three subprograms 
are GYCPS which would Initiate both files, GYCRP which would 
write to the running file, and GYCPP which would write to the 
profile file. 
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7. INPUT AND OUTPUT 


Table 1 lists the variables read on each input card along 
with the format for each. Definitions of the variables are 
included in the Glossary (Appendix P). 

There must be at least two and no more than two hundred 
Cards l6.n (see Table 2). They produce an Internal table 
specifying the velocity, enthalpy, and pressure at the edge of 
the boundary layer and the wall blowing velocity and either 
the enthalpy or enthalpy gradient (heat transfer) at the wall, 
all as functions of downstream distance x . The program uses 
linear interpolation in this table. For values of x outside 
the range- of the table, it extrapolates using the first two or 
last two entries. The input values of x must be in ascending 
order. The end of the table is flagged by a nonincreasing 
value for x or by n reaching 200. In the former case, the 
values on the terminating card are not used. If all three 
quantities are constant with x , only two cards, including 
the terminating card, need be used; the program will automati- 
cally generate another point to store in its table for inter- 
polation. As noted at the foot of Table 2, hg and^ Pg need 
not be input. (An exception occurs for pg when GMOMS = 0 . 

In this case, the values supplied for Pg will be used whether 
positive, zero, or negative.) When the program needs them at 
a particular x , the interpolated values from the Internal 
table are first determined. If these are zero (or negative), 
value, s computed from the homentropic relations for a perfect gas 
are substituted. (If the inputs HSTAG = hg and PSTAG = Pg 
are not needed for these calculations, they are not used at all.) 

There must be at least five and no more than JMAX Cards 
31. j (see Table 3). The values of y must be in ascending 
order. The end of the set is flagged by a nonincreasing value 
of y or by j reaching JMAX. In the former case, the values 
on the terminating card are not used. 
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Provision is made for specifying the initial profile of A 
on Cards 31* j- This is intended chiefly for use when A is 
solved for by a differential equation in the future. If zero 
is supplied for vA ^ the values will be determined in the usual 
way from the scale parameters (see Note 1 of the Glossary, 

Appendix F). . 

Cards 32 and 33 are not read at all if the flag terminating 
the set of Cards 31 -j is negative. The corresponding dependent 
variables (see Table 4) are initialized to zero for all values 
of y in that case. If the flag terminating Cards 31-0 is zero 
or positive, then Cards 32 and 33 are read. The end of the 
set is flagged by the first value of Y scaled by GV(1) which 
fails to equal the corresponding value of Y scaled by FV(1) 
in Cards 31. jj taken in order, or by J reaching the total 
number of cards in that set. In the former case, the corre- 
sponding dependent variables are initialized to zero for the 
remaining values of y . 

As' Indicated in Tables 2, 3 , and 4, the card sets designated 
l6.n, 31. j, and 33 -j are multiplied, or scaled, by elements of 
the vectors EV, FV, and GV, respectively. It is recommended 
that the beginner set all values of these vectors equal to one 
and forget them until a need is felt for them. Two typical 
situations where they would be useful are: (1) The Initial con- 
ditions are supplied as a function of y in, say, inches, 
and it is desired to use feet in the program. Then PV(1) (and 
GV(1) if needed) could be set at .0833333 and y values in 
inches punched in Cards 31. j (and 33. j if needed). (2) After 
one run starting with a certain level of turbulence, it is 
decided to make another under the same conditions except that 
the initial turbulence level is different, or zero. Then it is 
only necessary to change one card to prepare the input to the 
new run. 
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Input may be in. any consistent system of units; outputs 
will be the same units. Nondimensional inputs may be used, 
producing nondimensional output ^ but they must be consistent. 

For example, if BLAMO is made nondimensional with respect to 
a reference length, , then all other input lengths (DELTI, 

DTMAX, etc.) must be supplied as ratios to and all output 

lengths (X, Y, R, etc.) will be referred to . The intro- 

duction to the Glossary contains additional information on 
dimensions and nondimensionalization. 

A deck for one job might be set up as follows: 

Card 1 with INPLG = 1 

Cards 2-31. j, input for first run 

Card 1 with INPLG = 2 

Cards 2-29, restart input for first run 
Card 1 with INPLG = 1 
Cards 2-31. j, input for second run 
Card 1 with INPLG = 0 

This assumes a negative flag terminating the profile input is 
included in Cards 31- j- The above job, then, consists of two 
completely separate runs, the first being restarted after its 
first halt. The restart could be used to introduce a different 
cone angle, for example. If a run which halts on an error 
condition is set up for a restart, the restart probably won't 
be what’s intended but the error condition is likely to recur 
soon. On the other hand, it's possible to design the changes 
introduced on restart precisely to overcome an anticipated error 
condition. See Note 7 of the Glossary (Appendix F) for some 
specific information on restarts. 

Notes 5 and 6 of the Glossary contain information on inputs 
for special modes of operation. 
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Output consists of title pages, minor printouts, and major 
printouts. The title pages display the information input on 
Cards 1-29, for the most part in the same order. The labels 
are included in the Glossary, but one item near the top of the 
second title page is not labeled. Normally this is F but if 
the option of zeroing out correlations involving p’ and h’ 
has been selected, a T appears there. 

A minor printout gives basic information on the progress 
of the run and includes most of the auxiliary quantities described 
in Section 4. Two unlabeled items appear in the top line. On 
the left, there appears the name of the parabolic variable for 
which the smallest tentative Ax for the next step was computed, 
as described in Section 5* The letter Y appears there if such 
a calculation has not just been made, e.g., at the start of a 
run. On the right, there appears the name of a variable v/hich 
did not satisfy its edge tolerance until a point was added, as 
described in Section 3* It is not necessarily the only such 
variable. If no points have been added at the outer edge since 
the last minor printout, the letter Y appears there. 

A major printout consists of a minor printout plus profiles 
of the dependent variables and of the auxiliary quantities which 
are functions of y . 

A minor printout occurs alone only if NSTEP is a multiple 
of NIOLP. A major printout occurs: (1) if NSTEP is a multiple 
of NPOLP, (2) at time breaks (TBRKV), (3) if the maximum or the 
percent spread has changed by more than PCPMX or PCPSP times 
the value at the last major printout, (4) if NSTEP reaches 
NSTMX, or (5) if an error condition occurs. If NPOLP is one 
Oa major printout at every step is specified), two extras are. 
included: (1) profiles are printed after points have been 
inserted in or deleted from the y mesh by GYCAA and (2) a 
major printout shows the discarded results generated before 
each backup. 



A run halts when the. absolute value of a negative time 
break is reached, when NSTEP reaches NSTMX, or when an error 
condition occurs. 

Error conditions produce messages indicating the nature 
of the problem. A.R.A.P. should be consulted if a cure doesn't 
suggest itself. 



8. PINAL REMARKS 


Along with the GYC progranij an "SSP" listing of it was 
delivered to ARC. The SSP listing (named after the A.R.A.P. 
program that produced it) is an alphabetical index of variable 
names and subprogram names, locating every reference to each 
by subprogram and position within it, and printing the lines 
so located. Since GYC makes much use of COMMON, the SSP 
listing is a valuable aid in tracing the modification and 
uses of variables. 

Anyone who wants to understand the details of the oper- 
ation of GYC must turn to the program listings. Considerable 
effort has been expended (with varying degrees of success) to 
reduce clutter and to make the listings readable. With rare 
exceptions, the program flow is from top to bottom. Statement 
numbers appear in order according to an obvious -scheme . DO 
loops are indented and occasionally indentation is used to 
set off other statements. Natural breaks in the program are 
signalled by blank lines. (Rows of asterisks or other symbols 
are not used. It has been realized for several hundred years 
that decorated pages are not easy to read, attractive though 
they may be . ) 

With the descriptions of Section 6 as a start and with 
the aid of the SSP listing, those familiar with Fortran should 
be able to determine what is going on throughout GYC. 
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APPENDIX A.- THE EQUATIONS BEFORE MODELING 
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= -u ! p ’ . - u I p ’ . + 
1 >J J ,1 


" ("j!^ " "^!o) ^ " ¥H) * 


, ]ujy’ + 


+ .p«'ul + ,P*'U' + U*lu!u';‘. + u:u’“‘ j + 


,m 


,m 


+ y^.ulu'”^ + u^.u'u'^ 
j ,m ,j 1 ,m 
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+(ph'u'*uj + (p'u'*']lh'uj + 


+ + p'h'u\^j + p(u'Vh^^ + h'u''‘u.^J + 


+ p'u'^ulh p + p’h’u'^u. . 


= u'p^ + + p.jp'M + p'"pip;i - h'p;i + 


£ 


(i 


£ . — £ 


+ (MSr + PC~„ + p'u:>*' + p'u'^)u^u’ 


m 


,m m 


m 

£ 


+ + v*’u;j; ) u’u' ^ + Suyuiu-” + u!u;“ ) + 


+ v(“r + “!iiv >^*'4 ^ 

- ic „T'ul’^ + T»!'k'u! + T + 


,£ i 


,£ i ,£ 1' 


+ s(h'u”J P h'u'j> 5,,(h’u’>‘‘ + h’u;J) + 


+ ( 4 ,H + 4 ,i)>''p”*' ( 4 ,ii 4 , 1 )’“ 




m 
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pu^^h’h' ^ ^ + ^pu'^h’h' + u^p'h’h’) ^ - ^p'u'^^^h’h.' + 


"rrS-r 


+ 2p'h’u’^h » + 2ph .u'^h’ + 2h p'u’^h’ = 




yV 


I 


= 2h'pl + 2u h'p'j^ + 2p ^h’u’ + 2h'u’ p% + 


I 


. , ^ , T , I \ 1 - I J. 

+ 2fuu^ +yu^j^ + y’u- + U j h-u^^ + 


m 


+ + y*'u;{J) h'u'^ + sssy^'u””' f h'u;™) + 


y*'h' + 2y^u^„h»u»™ 
,in jX, ,m 


+ VI 'h' + 2u~„u"'^y*'h' + 2y^u~oh»u’;: + 


+ 2kh’T'J^ + 2k ^h*T”^ + 2T’Jk’h' + 2T^^h'k'’’ 




-£- 


-Z~ 


I 


, 5 - , — 

' 11 ^ r\ 


+ p'p' u + p'uj^ u + p’u’ ■'■ 


+ P'P'vi'^Uj_^^ + p'u'^uj^ p^^ 


- -p'p^. - pSu^u'^ - 2B(u^j,p'ui + P'uiu;j) - >^‘iP'P'u: + 


+ p(p'u”J + P'u;Jj) + 5^i(p’4’*' + p'«;J) + 


♦ {h,i. * “t,i)p'p”‘^ + (“i,ii 


+ p'y*i + u“ -y*'p' + y*p'u'^'‘. + y^.p'u’ 


-m 


Jti 


,m 


,m ,1 


,mi 


,mi ,m 
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APPENDIX B.- THE MODELS 


uj^ujUjJ. = -VUUqA^ 

j 

(M“j ) 

t 

ujujp' = -VRUqAf 

/ 

(p'u') 

k 

u!u!h' = -VUHqAl 
1 J \ 


u^pTpt = -VRRqA 

P’P' j 

u|p'h' = -VRHqA 

Q. 

u^h’h' = -VHHqA 

(^h’h’ ) 


ul 

1 J 


I I 

u! ’ h' 

1 ,m 


^u'u! + b eulul + (l-6)g . — ~ 
^ ^ J A,. ^ J -^J ^ 


I ^ t 

u' u' 


A- - J 
^AHU 


AV 


+ BHU ^ h • u ! 
Ay ' ^ 


P m 
^ ^ m 


h' 


,m 


ARU 

a 2 

AHH 

a 2 


+ BRU 


£q 


f11 ' 


Ay 

+ BHH ^ 
Ay 


p 'u 


h'h' 


, ,Z , 


P £ 

where it„ = -BDV — ^ 

YP 


Z 



p ' u’ = ir p p ' u' 


p’u^Iu’J = VBSit^P’u|u’ 




,k 


“Ui 


, k .1 

TT, 7T. U' U’ 

k £ 




P ^ 2^l[^ ^ ),± 

^ = K^li 


ulul 


i J,k 




2 (®Jk"i'*U ®lk’^J“ln) 


u! .u': = WMD u 

1>3 V 


k;i),d - i®ij("’‘'";'),k " i8i3kk“: 


g .u!u' 
^mj 1 n 


^ g ulul + 
3 ®nm 1 j 
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= (wWGUq/A3)(uiu- - i 


, ,m ,n 
u!u' u’ 

1 ,n ,m 


= (q/A) WVJUl(|u^u'^) + WVIU2 uju'^ 


,m , n ,S, 
u* u' u'n 
,m 


= (q/A) WWDl(u’\’^) ^ + WWD2 u’j^u'J 


h'u'%'^ = (q/A) 
,n ,m 


WWHl^ 


h'u’^) ^ + WWH2 h'u’J 


piut^ ,n ^ MRlfp'u'^) . + WWR2 p’u’J 


p'u! = pA^ PMau"^ ulu'^ + PMU2 ulu'%’^ 
1 ,n 1 ,m 1 ,n ,m 

p*h' = pA^ PMHu”^ h'u'!^ + PMH2 h’u’%’^ 

,n ,m ,n jin 


P' P’ 


pA^ PMRu^ p'u'’^ + PMR2 p’u’%'’^' 

^ jH jH jin 




pA^ PGUu"^ u! .u’^ + PGU2 ul 

' ,n ,m i,j ,n ,m 


p'h’, = PGHp WGHlu“ h'ul + WGH2u“‘ h'u’ + WGH3u. „ h'u' 
,1 ,m 1 m i,m 


+ PGH2«WWGH ^ h’u^ 


p'p'. = PGRp WGRlu ‘ p'u' + wGR2u . p’u' + WGR3u. ^ p'u*' 

jiiil jj-ni -L^iJi 


+ PGR2^^WWGR ^ p'u! 

A 1 
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U'p' + uVp>^ = PTUM ((SJp^)^ + 5 ^(u!_p’)^5,) 

h'p^ + u^h'p'^ = PTHM ((h'p')^ + ^^(h'p')^j,) 

p'h' = -h’h' p/h 

p’p* = -p’h’ p/h 

p ' h ’ h ' = 0 

u ' ! p ’ = 0 
^ ’Jl 

h'u'^p'^ = 0 

The following constraints apply: 

3-WGHl + WGH2 + WGH3 = 0 
3*WGR1 + WGR2 + WGR3 = 0 . 
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APPENDIX C.- THE EQUATIONS AFTER MODELING 


» 


Mean continuity 


(RS-vUS"L) »L + RU"LrL = Q 


4 


mean momentum 


RS+US"L*US« I ,L + RU"L»US'I,L 
+ {.RS*UUt'L*I + US"L*RU»I),L 
- VRU*(0L*(RU* I IL + RU"LtI))»L 
= _ PS, I + (EMU+LIS'IIL + EMU*US"L,I ) ,L 
+ EMUT=t:WMT+{ rU* I IL + TU"L,I),L 
~ 2/3*EMUT* ( WMT*TU"MtM - BDGP^TU "M^PS , M ) , I 
+ <EMS*u&M,v,,fq f EMST*B0GP*1U''M*PS,M1 , 1 


4 


MEAN enthalpy 


RS*US"L*hS,L + RU"L+HS,L + <RS^HU”L + US"L*RH),L 
- VRH*(QL=»-RH,L) I L 
= bMOMS*( US"L*PS,L 

+ ( RS + BLAM + BI. AM+- ( PMU + US”M,M*( 1/2*UU"LN,M 

- l/4*(ai«N'M*UU"LK,K - S "L ’ M^UU"NK , K ) 

+ 1/2*(&"N’M*BDGP*UU"LK*PS,K •* i"L*M*BOGP + UU“NK*PS»K) ) 

+ PMU2*0/BLAM*(WWU1>KUU"KL,K + WWU2*BDGP=»=UUiiKL*PS,K) )),L 
- RS*8LAM+BLAM*( PGU*US"M , N* { UMD* ( ( BDGP*UU" NL*PS » L ) ♦ M 

- l/3*&"N*M*(BDGP*UU"KLi'PS,L) »K ) 

+ 1/3*S"iM'M*BDGP*BUGP*UU"KL*PS,K*PS,L) 

+ PGU2*Q/BLAM*(WWD1*(BDGP*UU"KL*PS,K) ,L 
+ WWD2*bDGP*B0Gp)KUU"KL + PS,K*pS,L) ) 

+ (EMU*{US*M!L + US"L,M) 

+ EMUT*( WMT*{TU*M1L + TU"L,M - 2/3*a''L » M*=TU"K » K ) 

+ 2/3*&«L»M*B0GP*TU"K*PS,K ))*US"M,L 
+ (EMS*US"K,K + EMST*BDGP*TU"K*PS,K)*US"L,L 
+ EMU*A/BLAM/3LAM>i=UU''L*L + B/BLAM*RS*Q*UU"L*L 
+ EMU*(UU"LM,LM - 2*'< BDGP*UU"ML*PS , L ) , M 
+ BOGP*BUGP:i^UU"ML*PS,M*PS,L) 

+ emUT*(US'L!M + US"f 1 »L) + (WMT=t<( TU"L,M - 1 / 3 *S"L ’ M* TU"N , N ) 

+ l/3*a"L*M+BDGP*TU"N*PS,M) 

+ rMS*BOGP=t^BDSP>fUU'»MLTpS,M*PS,L 
+ emST+US"i'1,M*BOGP*TU"L»PS,L ) 

+ tCAP*HS,U!L + 1/2*CAYT*TT,L!-L 


oSf PAGE IS 
OOE QUAlttj? 
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uu 


RS*US''L*UU» IJ»L - ( VUU*RQL* (UU* IJ ! L + UU"L*JiI + UU"L»I,J) 

+ VRU*QL*US"L%( RU’ Ii J + RU*JtI)),L 

- RU*'L f L*UU* IJ + RU* I + US"L*US» JtL + RU • J*US'< L*US * I ♦ L 
+ RSfc(UU''L»l*US' J«L + UU"L’ J+US'i ,L) 

- VRlj*OL»{US» J,l.*:(RU* I IL + RU"L,I) + US * I , L% ( RU » J! L + RU"L»J)) 

= -( RS+BLAnH.8LAM*(F’P1U*US"MiN*( 1/ 2*UU»M » I . M 

- - n1» TM*UU"LIM»L) 

+ 1/2*{&''w«M*BDGP*UU"L* I*PS*L - 51 ' IM*BDGP*UU"LN4:PS » L ) ) 

+ PMU2*Q/BLAM'*(WWUl*UUiiL»I,L + WWU2*B0GP*'UU'*L » I*PS , L ) ) )»J 

- ( RS*BLAM+BLAM* (PMU + US"M*i\i* ( 1/2 + lJU ”M* J ♦ M 

- l/‘t*(i"N'M*UU"L*JiL - a* JMi|cUU"LNil_) 

+ l/2*(gMM«H*8DGP*UU"L» J+PS,L - 3 ’ JM*BDGP*UU"LN*PS » L) ) 

+ PMU2+Q/BLArF*lUWUl*UU"L' J,L + WWU2*B0GP»UU”L • J*PS , L ) ) )»I 

+ RS+( PGU*WGU+=US"M»N*(a’MJ*UU"l\i» I + 3 » M I *UU"M ♦ J 

- 2/3*S"N»M*UU» IJ - 2/3 + a * I M + 2/9*3 * IJ*g'' M • M*UU»L ♦ L ) 

- 0/BLAM*(UU*IJ - l/3*a»IJ*UU"L*L) ) 

+ (FMU*UU*XJ,U lU 

- 2+EMU*A/BLAM/BLAM*UU* I J 

- c>*B/BLAM*RS*a*{BETA*UU» IJ + 0M8ET*a • U/ 3 *UU"L ♦ L ) 

+ IEMU + EMS)*( (BDGP*UU"L* I*PS,L) , J + ( BOGP*UU»L * J*PS ♦ L ) , i 

- WMD*( (B0GP*UU"L'I*PS,L) , J + ( BDGP*UU"L * J+PS , L) i I 

- 2/3*a*IJ*(BDGP*UU"LK*rS,L) »K) 

+ 2/3*a» IJ*BDGP*-8DGP*UU''LK*PStL*PS,K ) 

+ rMU,L*( 1/2*UU"L*I»J - 1/4*(&"L» J*UU"N* I .M - 3 » I J*UU»LM , M ) 

+ 1/2*(&'’L* J:*BDGP*UU"H»I*PS«I'1 - a ' I J*RDGP*UU»L''1*PS , M ) 

+ 1/2*UU"L'J.I - 1/4*(S>(L' J,M . 3'IJ*UU"LM*M) 

+ 1 /2*(S"L'1*80GP+UU"M» J*PS,M ~ 3 ' I J*BBGP*UU"LM*PS , M) ) 

+ q:MUT*{ {JS*I,L + US» L, I ) *(TU* J! L 

- WMT*{TU»JiL - X/3*2,"L» J*TU"niM) 

- 1/3*&"L* J*BDGP*TU('M*PS,M) 

+ (US*J,L + US«L,J)*(TU*l!L 

- Wf1T*(TU’l!L - 1/3*&"L» I*TU"M»M) 

- ,l/3+g"L» I*BUGP*TU'iM*PS»M) 

+ TU’J*(US»1»L!L + US ’Ll I !U + TU » I * ( US • J , L ! L + US’L,J!L) ) 

+ rnSl*( US"M,M*(TU» I,U + TU'J,i - WMT*(TU*I*J + TU*J,I 

- 2/3*a» I J*TU"L,L) + 2/3*a* Ij*BDGP*TU"L*PS,L) 

+ tU» J*US"M, ril + TU*l*US"ri,MJ ) 

+ ri«)SiI*8DGP*UU"’1‘ J*PSiH + EMS ♦ J*80GP*UU"M • I*PS t M 


SffSSS 
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ft 


HU 


RS*US''L*HU» I ,L - (VUH*RQL*(HU*I !L + HU''L»I) + VRH*QL*US"L*RH , I ) , L 

- RU"L«L*HU*X + RU'I + US"L*HSiL + RH + US " L*US * I . L 
+ RS* (UU"L* I*HS,L + US'-I *L*HU"'U 

- \/RU*QL*(HU'X!L + RU"L.I)*HSfL - VRH:«cQL*RH ! L*US * 1 1 L 

= SMO?lS+( PTUH*US"L+ ( RS + BL AM+BL A * ( P-'^U*US"M « N + l/a+UUi'N * I » M 

- 1/H*PMU*(US'»1^,M + UU‘'K» I,K - US » 1 1 W*UU"NK , K ) 

+ 1/2*PHU*(US"M,M*30GP*UU"K' I*PSiK 

- US*X»N+BDGP*UU"K[\i + PS»K) ) 

+ PMU2*KQL*(WWU1+UU"K» I iK + WWU2+BDGP*UU"K * I*PS , K > )»L 
+ UU"L'I*Pb«L ) 

- ( RS*BLAM*BLAM* (PMH*US"M»H+ { WMH=«c ( HU " M i FI 

- l/3*i"N»Fl + HU"LtL) + 1/3=^&"W'F1*30 GP*HU"L*PSiL ) 

+ PF1H2 + 0/BLAM* ( WWH1*HU"L «L + WUH2*30GP*HU" L*PS » L ) ) )fl 
+ P6H*RS*( WGHl*US"Fl,Fi*HU» I + UGH?*US " 1 I + HU • M 
+ WGH3+US»l*Fl*HU'iFi) + PGH2+WWGH*KS*Q/BLAi«l*HU » I 
+ GFI0MS*( (EMU^{US»F1! L + US"l,F1) + EMUT/CP*( WFiT* { HU ♦ FI ! L 
+ HU"LtM - 2/3 *&"L'H*hU"K,K) 

+ 2/3*S"L’F1=^BDGP*HU"K*PS*K ) ) + ( 1/2*UU"M’ I » L 

- l/4+(&"Fl*L*UU"N>I ,N - a'lL*UU"MN»N) 

+ l/2*{S''Fi*L*B0GP*UU''N»I*PSiN - 3 • IL*BOGP*UU*'FllM*PS ♦ N ) ) 

+ ( EFiS+US"K iK + EF1ST/CP*BDGP*HU'»K*PS 1 K ) *BDGP*UU"L* I*PStL 
+ EMU+US"LiM*( 1/2*(UU* XLIFI + UU"Ft'l,L) 

- 1/4+ (2+a,"Fi’L + UU"K» I ,K - S"Fi • I + UU"K * L t K - 3 * XL-*UU”KF » K ) 

+ 1/2* (2*&"F1'L*BDGP*UU"K * I*PS tK - g,"M * I *BDGP*UU"K ’ L*PS t K 

- 3*IL + BDGp*UU"KFl*PS,K) ) 

+ EF1UT/CP*US"L,M*(US»LIF1 + US"F!fU+HU’I 

+ EF1ST/CP*US"L,L*US'‘MiF 1*HU» I + EF1S*US"L» L*8DGP*UU»M* I*PSi V) ) 

+ (CAP*HU»I,L) !L - CAP*hUDFlF*HU* t 

- {CAP*(WHT*(HU*I !L - 1/3*&"L* I*HU"K,K) 

+ 1/3*X"L» X*aDGP*HU'tK*PS »K ) ) 1 L 

+ CAYT/CP/CP!<^{ {hU*I*HStU !L - HS.L*( v^FlT* ( HU • I ! L 

- 1/3*&"L'I*HU"K,K) + l/3*X"L' I*BOGP*HU»'K*PStR )) 

+ EFiU*(WMH+(HU*1iL - l/3*3’lL*HU"KiK) 

+ 1/3*3* IL*BdGP*HU"K*PS,K) !L 

- EPU*HUUF1F*HU* 1 + 1/?*IEF1U + EFIS ) * ( 3DGP*HU»K*PS i K ) » I 
+ EFUtL*( WF1H*(HU* I I L + HU"L»I - ?/3 + 2." L ' X *HU "K t K ) 

+ ?/3*&"L* l*BOGP*HU"K*PStK) 

+ EF.UT/CP*( (US*I,L + US’L,I)*1/2*HH1L 
+ (US*I,L + US*U,I) 1L*HH) 

+ EF.ST/CP*(US"M,.'1*1/2*HH,I + US"F1,MI*HH) 

+ EFlSt I*BDGP*HU"F1*PS*M 


0® ?00B 
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HH 


RS*US"L*HHiL - VHH=t=(RQL*HH»L) ! L - RU"LiL*HH 

+ 2*RH*US»L+HSiL + 2 + RSi^HS , L*HU"L - 2*\/RH*QL*HS , L*RH 1 L 
= 2+Gt’10MS*{ PTh^l+US"L* ( RR*8LAM*BLAM* ( PMh*US"i'' 1 t ( WMH* 
(HU"N,M - 1/3*&"N'M*HU"K»K) 

+ 1/3*&"N' M*BDGP*HU"K*PS, K ) 

+ PMH2+Q/BLAM* { WWH1*HU"K»K + WWH2*BBGP*HU"K*PS t K )) )tU 
+ PS*L*HU"L 

+ (Enu=t:( US’MIL + US"L,M ) 

+ EWUT*( WHT*(rU'i^!L + TU"L,''^ - 2/3*&"L ' M*TU"K» K ) 

+ 2/3=t=^"L' M*BDGP*TU"K + PS 1 K riU'iMiU 

- l/3*ii"M*L=»-HU"W» N ) + l/3+X"M * L*BOGP*HU"N*PS i N ) 

+ {EHS*US"K*K + Ei^ST+3DGP*TU"K:t:PStK)*B0GP+HU"L*PS,L 
+ EnU*US"LtM*(WMH*f HU'LIM + HLM'MiL - 2/3*&"M ♦ L*HU"K , K ) 

+ 2/3 + &"M'L*BDGP>I'HUi'K#PSiK) 

+ EMUT/CP*US"L,M*(US*L!M + US"M.L)*HH 

+ eMST/CP+US"L*L*US"M*M*HH + EMS+US"L»L+BDGP*HU"M*PS»M ) 

+ (CAP*HH,L)!L ” 2*CAP*HHDi^F*HH 
+ CAYT/CP/CP*(2=tHS,L!L*HH + HS,L*HH1U 


u RU 


RS*US"L+RU’ I,L - (VRU*R 0 L*(RU»I!L + RU"L,I) 

+ \ZRR*QL*US"L + RRi 1 ) ,L - RU"L.L*RU*I 
+ RR+US"L*US’ 1 tL + US"L*RS.L*RU’I 
+ RS*US’I,L*RU"L + RS*RSiL+UU"L * 1 

- \/RR*QL*RR!L + US*ItL - VRU*QL*RS , L* ( RU » I ! L + RU»Lil) 

= -CRS*BLAM + BLAH* (PFiR*(WHR*(US"Mff\)*RU"I'JiM - 1 /3+US "M , i^*KLJ " L . L ) 
+ 1/3*US"M,M*BDGP*RU"L*PS,L) 

+ PMR 2 *Q/BLAM=+ (WWRX*RU"LtL + WWR 2 : 4 ^ 3 DGP*RU"L*PS i L ) ) ) » I 
+ PrR*RS+(WGR1+US"M,M*RU*I + WGR2*US"MiI*RU'M 
+ WGR 3 *US*iiM*RU"ri) + PGR 2 *WWGR + RS + Q/BLA' 1 *RJ» 1 
- RS+RS*BDGP*UU"L' I*PStL - 2 *RS+US"L»L*RU* I 
+ 2 *VBS*VRU*R 0 L+ 3 UGP+pStL*(RU'I 1 L + RU"LiI) 

+ VRR*QL*US"L,L*RRi I 

+ EFU*(WMR*(RU’ IiL - 1 / 3 * 3 ' lL*RU"KfK) 

+ l/ 3 *a* IL*BDGP*RU»K*PS,K) !L - EMU*RUDMF*RU’ I 
+ 1 / 2 *(EMU + EnS)*(BDGP*RU"L+PStU »I 
+ Er!UtL*(Wi^R*{RU'I!L + RU"L,I- 2 / 3 *i" L ' I *RU"K , K ) 

+ 2/3*&"L ' I*BDGP*RU"K*PS, K) 

+ EMUT*( (US' I ,L + US'Ltl)*l/ 2 *RTIL + (US'I.L + US'L»I)lL*Rn 
+ EF;ST*(US"M«M* 1 / 2 *RT,I + USi'M,MI*RT) 

+ El^Si I*BDGP*RU"f'^*PSiM 


OEIGWAI-Pf-®® 
OE POOR QUAIiTES 
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APPENDIX D.- THE PINAL EQUATIONS 


fl FIRST PASS 


RSUS*(R+WS) :i 
+ RSVS* ( R'*WS) ; 2 
+ Rv<-(R*WS):2 
+ (rs*R*vw ):2 
=: ( fMU* ( R*WS) ; 2 ) ; 2 


RSUS*(R*VW) ;i 
+ RSVS+(R*VW) ;2 

- 2 *r+RS*Ww*wS 

- 2*(VUU*RQL*(R+y/W) :2) ;2 

- RV2*R*VW 

+ US*RV*(H*Wb) J1 
+ VS*RV*(f<*WS) :2 

- C*iaIS*RW + WS 

+ Rs*uv* (R*ws ) : 1 
+ Rs+vv* ( R*ws ) : 2 
+ RS*R*VW*VS2 

- 2*\/RU*QL*R\/2+( p*kiS) ;2 

= /!4*PNU*R=^R*\/V2*(RBLS*(1/R*WS) ;2) ;2 

- 1 /14.*PMu*RBLS+WSR2*(R*H*{ VV) :2) :2 

- fPMU2*WWUl*RQL*(R*VW) :2) ;2 

+ pgU*WGU*RS*R<^R*UV*( 1/R*WS) ;i 
+ PGU*WGU*RS*R*R*\/V* ( 1/R*ws ) : 2 
+ 1 /3^S*PGU*WGU*RSUS*Vht 
+ S + pGU*WGU*RS*U\/+WS 
+ c*PGU*WGU+RS*VV*WS 
+ 1 /3»PGU*WGU=»-RS*\/S2*K*VW 

- c*pGU*lrJGU*RS*WS*WW 

- 5/3*PGU*WGU*RS*US1*R*VW 

RQOL*K*Vs«i 

+ (emu*(r*vW) : 2) :a 

- 2^A + L'MU/8LaM/BLAM*K + VW 

- 9*8*BETA>t:RQCL*R*VW 

+ 1 /A+EMua+tRifVW) !2 

+ (1 - Wmt/5)*EMUT+TV2*(R+WS) ;2 

+ emut*tv*( (R+ ws) :2) ;2 


OB BOOB QOAUTY 
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)4 


SECOMD Pi\SS 


RSUS*(US) ;i 

+ RSVS*(US);2 
» S*RS/R*1/JS>»:vJS 
+ RV*(US) ;2 

+ (RS*uv):2 

= -nPDX 

+ ( emu* (US) ; 2 ) : 2 


Rsus*(uv) ;i 

+ RSVS*(UV):2 

- 2*S*RS/R*WS*VW 

- 2*r*RS/R*WS*UW 

- 2*(Vuu*rql*(uv) ;2) :2 

- RV?*UV 

- C/R*WS*WS*RU 
+ us*R\/* (US) : 1 
+ VS*KV*(US);2 

- S/R*WS*WS*RV 
4 RS*VS2*UV 

4- RS*USl*Uv 
4 - RS*VV*(US) ;2 

- 2*\/RU*QL=I‘R\/2*( US) :2 

= -1 /4*PMU*V\/2*(RRLS*(US) ;2) :2 

- t/4*PMU*R3LS*US2*( (VV) :?.) :2 

- (PMU2*WWU1*RQL*(UV) :2) :2 
4- i/3*PGU*iJGU*RS*VS2*UV 

- c*PGU*WGU+RS/R*WS*UW 
+ 1 /5*pGU*WGU*RS*USl*U\/ 

4- PGU*WGU*RS*VV*(US) :2 

- s+pGU*WGU*RS/R*WS*\/W 

- 2/S*S*P6U*WGU*RS/R,*US*UV 

- ROOL*U\/ 

+ (EMU*(UV) ;2) ;2 

- p+A+EMU/BLAM/BLAMcUV 

- ?*R*BETA + RQ0L*'J\/ 

4- 1 /(4*EMU2+ ( UV.) : 2 

+ (1 - WMT/3) *EMUT*TV2* ( US) ; 2 
4- rML)T*TV*{ (US) :2) :2 


RSus*(R*uW) :i 

4 RSVS*(R*UW) :2 

- 2*S*RS*WS*wW 

- ( vuu*RQL*(r*uw) :2) :2 

- RVP*R*UW 
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+ US*RU%RWS1 
+ VS*-RU*RWS2 
+ R*IJS*RW* ( US) : 1 
+ R*\/S*RW*(US) :2 

- S*WS+WS*RW 
+ RS*UU*HWS1 
+ RS*UV*RWS2 

+ RS*R*US1*UW 
+ RS*R*VW*(US) 12 

- VRU*QL*RWS2*(PU) ;2 

- VRtJ*0L*US2*(R*RW) J2 

= PGU*WGU*RS*WSR1*R*R*UU 
+ PGU+WGU*HS*WSR2*R+K*UV 
+ 1 /3*S*PGU*WGU=t=RSUS*UW 
+ S=»-PGU*WGU*RS*-WS*UU 
+ c*PGU*WGU*RS»:WS*U\/ 

+ 1 /3*PGU*WGU*RS*USl+h*UW 
+ PGU*WGU*RS*U32*R*VW 

- S*PGU»WGU*-RS*WS*WW 

- P/3+PGU*a/GU*RS*\/S2*R*UW 

- rqol*R*uw 

+ (EMU+(R*UW) : 2) :2 

?*A*EMU/8LAM/BLAM*R*UW 

- p^R+bETA^^ROOI + R + UW 

+ <1 ■* WMT)*Ef^UT*RTW2*(US) ;2 
+ (1 - 1«MT)*EI«’UT*RW&2*1TU) ;2 
+ FnuT*R+Tvi*( ( US) J 2 ) :2 
+ n^UT*TU*RWS 22 


0 ® 
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H 


third PA.SS 


Rsus*(uu) :i 

+ rsws*(UU):2 

- 1|.4S=I'RS/K*WS*UW J 

- (V1JU*RQL*(UU) ;?) ;2 

- RV?*UU 

+ 2+US*RU*US1 
+ 2*VS*RU*US2 

- 2*S/R+WS+WS*RU 
+ 2*RR*HS1+UU 

+ 2*RS*US2*UV J 

- 2*VRU*QL*US2*(RU) :2 

= 8/<5*PGU*wgU»RS*US1*UU 
+ a/3+PGU*WGU*RS+US2*UVJ 

- p+S*PGU*WGU*RS/R*WS*UWJ 

- 4 /g*PGU*wGU*RS+VS 2 *UU 

- u./9*S*PGU*WGU*RS/R + US*UU 

- 4/9*PGU*W6U*RS*VS2*VV 

- ?/3*P<5U+iJGU*RS*WSR1*R*UWJ 

- p/3*PGU*WGUN'RS*WSR2*R*VWJ 

- 4/9+S^PGU*WGU*RS/R*US*WVi 
+■ p/g + PGU*WGU+RS*USl=*=VV 

+ p/9 + pGU*vJGU*RS+USl*WW 
+ ?/9*PGU*WGU*RS*'VS2*WW 
+ 2/9*S*PGU*WGU*RS/R*US*VV 

- o/3*RG0L+UU 
+ 1/3*RG0L*\/V 
+ 1/3*RS0L*WW 

+ (EfTU4--(UU) ;2) ;2 

- 9*A*EMU/3LAF1/BLA(’UUU 

- 2»b*be;ta+rqol*'Ju 

- 2/3*B*0MBET*RQ0L*UU 

- ?/3*B»OM3ET:fRQ0L + \/V 

- p/3*8*0MB£T*RQ0L*t;W 

+ 1 /2*EHU2*{ VV) :2 

+ ?*(1 - W(’1T)*E:MUT*US2*{TU) ;2 
+ 9*EMUT*TU*US22, 


RSuS* ( yv ) : 1 

+ RSVS*{ VV) t2 

- 4 + f:*RS/R*vJS*VW J 

- 3*(\/UU+RQL*< vv) : 2) :2 

- RVp+VV 

- ?*r./R*WS*WS*RV 
+ 2 *RR*Vs2*VV 

= -P!«=(PMU2=»:WWU1<‘RQL*(VV),:2) ;2 
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+ tx/q*pSU*(fl/GU^RS+VS2*VV 

- ?5(.C*P<5U*WGU*RS/R*WS*VWJ 

- ii/g*PGU*WGU*RS^USl*\/V 

- 4/9*S=*'PGU*WGU*RS/R<'US*\/V 

- '4/9*pGU>t=wGU’^RS^LJSl*UU 

- ?/3*PGU*wGU*RS'*'USa*UVJ 

- 9/3*pGU*wiGU*RStWSRl*R*Ul-'J 

- :?/3*PGU*WGU*RS*WSR2*R*VWJ 

- 4 /g*S*PGU<=wGU*RS/R+US*W!(v 
+ ?/q»P&U*wGU*RS*USl+WW 

+ 9/9*PGU*wSU*RS^VS2*UU 
+ ?/9*PGU*wGU*RS*VS2^WUi 
+ ?/9*S*PGU'»^wGU*RS/R*US*UU 

- 5>/3*R00L + VV 
+ 1 /3*R0UL*UU 
+ 1 /3*R00L*I(JW 

+ vv) ; 2 ) :2 

- 2*A*E'MUy8LAM/BLAM+V\/ 

- ?*b*r£ta+RQol*\/v 

- ?/3*B*OM3ET*RGOL*UU 

- ?/3*B=*^0MBET*RQ0L*VV 

- 9/3*B*0MaET+RQ0L*WW 

+ F-Mua+ivv) :2 


F^SUS=K ( R*R*WW ) ; 1 

+ RSVS+(R*R+WW) ;2 

- (VUU*RQL.*(R*R*VJW) ;2) ;2 

- RV2*R*R*WW 

+ 2 f R*US+Rv/*RWSl 
+ 2*'R*\/S*F^W*RWS2 
+ 2)i^RS+R=*'RirfSi*UWU 
+ 2*«S*R*RWS2*\/WU 
-• 2*\/RU*QL*RwS2* ( R*RW I ; 2 

= fj/3*PGU*WGLIiiRS*WSRl*R*R*R+UWJ 
+ 4/3*PGU*WSU=*-RS*WSR2*R*R*R*VWJ 
+ r/9*S*P6U*WGU*RS*R*US*WW 
+ p+S*PGU*rtGU*RS*R*WS+UWJ 
+ p*C*PGU*t»GU*RS*R*WS*VWJ 

- 4/9*PGU*UGU*RG*USl*R*R*WU 

- 4/g*PGU*WGU*RS*VS2*R*R*WW 

- 4/g*pGU*WGU*RS*USl*R*R*UU 

- p/3*pGU*iJGU*RS*US2*R + R*UV 

- 4/g*pGU*tgGU*RS*VS2*R*R + VV 
+ 9/9 + PGU*v-GU*RS+US1*R*R*W 
+ p/9*PGU*vJGU=»:RS*\/S2*R*R*UU 
+ :?/9*S*^PGU*WGU*RS*R*US*UU 
+ :?/9*S*PGU*WGL)*RS*R*US*VV 

- o/3*RQ0L*R>kR*WW 
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+ 1 /3*R<30L*R*P*UU 

+ f emu* ( R*R*WW) ; 2) ;2 

- ?*A*ElMU/bLAP'i/BLAM*R*R*WW 

- ?*b*beita + rqol*r*r*ww 

- ?/3*B*0M3ET*RQ0L*R*K*UU 

- ?/3*B*0MBET*RQ0L+R*R*VV 

- 9/3*B*0MBET*RQ0L*R*R*WW 
+ 1/2*EMU2 + R*R*(\/V) ;2 

+ ?*(1 - WMT)*EMUT*RWS2*(R*TW) ;? 
+ 5*EMUT*R*TW*RWS22 
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4 fourth Pi^ss 

RSUb*(HS) :i 

+ RSVS*(HS) :2 
+ RV*(I1S) 52 
+ (KS*HV);2 
= &poMS*us*uPnx 

+ GMOMvS + eMU+US2*USJ2 
+ GMOMS*lMU*RWS2*WSRj2 
+ A*OMOMS*L^iU/BLAM/BLAM*OQJ 

+ b*omoms*hQol>^qqj 
+ ( CAP* (hs ) : 2 ) ; 2 


RSUS+CHV) :i 
+ RSVS*(HV);2 

- 2*C+RS/R*WS+HW 

- 2*(VUH*RQL*(HV) ;2) :2 

- 

+ US+RV* (HS) : 1 
+ \/S*F\V* (HS ) ; 2 
~ C/R*WS*WS*HH 
+ RS*\I\I*[HS) 12 
+ RS»VS2*HV 
- 2+VRU*QL+RV2*{HS) :2 
= GMOMS*UV+DPOX 

+ groms*v\/ + dpdy 

- (PMH2*WWHl*ROL*(HV) ;2) :2 
+ PGH*WGHi+rS+US1*HV 

+ PGH*WGH1+RS*VS2*HV 
+ S*P6H*WGHl*RS/R>t:US*H\/ 

+ PGH*WGH2+RS+US2+HU 
+ PGH*WGH2*HS*VS2*H\/ 

+ PGH*WGH2*RS*WSR2+R*HW 

+ C*P0H*iJGH2*RS/R*WS*hW 
+ PGH*WGH3*RS*VS2 + H\/ 

- C*PGH*WGH3+RS/R*WS*hW 
+ PGH2*WWGH*RQ0L*HV 

+ 3/2*GM0MS*LMU+US2*UvJ2 
+ 3/4+GP.0MS*£-MU=t.RWS2 + VWRJ2 
+ 3/A*GMOr'iS*EriU*WSR2=t=RVWJ2 
+ GrOi'1S*EMUT/CP*OS2*USJ2*H\/ 

+ GF-OMS+EMU [ /CP + WSR2*RWSU2*HV 
+ (1 - 2*Wi^T/3) +(CAP*(HV) :2) ;2 

- cap*hudmf=*=hv 

+ CAYT/CP/CP* (HV*(HS) :2) ;2 

- 2/3*CAYT/CP/CP*HS2*WMT*(HV) ;2 
+ 2/'3*EMU*(WMH+(HV) :2) :2 

- Er-U*HUDFF*HV 

+ I4/3*EMU2*WMH* ( HV ) ; 2 


rsus*(hh) :i 

+ RSVS*(HH) ;2 
~ vhh*(RQL*(HH) : 2) :2 
- RVf2*HH 
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+ 2*US*RH* (MS) :i 
+ 2+VS*RH+{HS) :2 
+ 2+RS*(HS) :2*HV 
- 2*VPH*QL*HS2*(RH) :2 
= 2*GM0MS*DP0^*HU 
+ a^GMOi’IS + DPUY + hV 
+ 4i^GMOnS*EMU*US2*WMH+(HU) :2 
+ 2tGM0MS=*-EMU*RWS2*WYiH* ( 1/R + HW) !2 
+ 24.G!^0MS*LMU*WSR2*WMH=»=(R + HW) :2 
+ 2*GM0MS+EMUT/CP*US2*US J2+HH 
+ 2=^G^^0HS*LMLJT/CP*HSR2*RWSJ2♦HH 

+ (cap*(hh) ;2) :2 

- 2*CAP*HHUHP*HH 
+ 2*CAYT/CP/CP*HH*( (HS) :2) :2 
+ CAYT/CP/CP*HS2*(HH) :2 


of PAGB is 

POOR QUALITY 
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FIFTH PASS 


RSuS+{t?U) :i 

+ RSVS*(HU ) 12 

- 2+S+RS/R*WS*RW 

- (Vru*ROL*(RU) :2) ;2 

- RV?*RU 

+ US+RR+USl 
+ VS*RR*US2 

- S/R*WS*WS*RR 
+ US 4 <RS 1 *HU 

+ VS=^r'S2*RU 
+ 3*RS*US1*RU 
+ RS*HS2*RV 

+ RS*RS 2 *U\/J 

- VRR*QL+RR2*US2 ' 

- VRiJ*QL*RS2*(RU) : ? 

= PfiR*WGKl*RS*USl*RU 
+ pGR*WGR1*RS*VS2*RU 
+ r*PGR*WGR 1*RS/R*US*RU 
+ PGR*WGR2+RS+US1*RU 
+ PGR*WGR2+RS*WSR1*R*KW 
+ S*PSR*WGR2+PS/R*WS*RIJ 
4- PGR*WGR34.RS*US]+RU 
+ PGR*VJGR3 + RS=tUS2*RV 

- S+PGR*WGR3*RS/R*WS*RW 
+ PGR2*WWGR*ROOL*=RU 

- P+RS+VS2+RU 

- 2*S*RS/R*US*RU 

+ emu*(Wmr* (RU) ;2i :.2 

- r|W!U+RUDMF*RU 

+ rMU2*WMR*(RU) : 2 
+ 1 /2*EMUT*US2*=( RT) :2 
+ EMUT*US22*RT 


RSus*(RV):i • 

+ RSVS*(RV);2 

- 24 r*RS/R+WS*RW 

- 2*{VRU+rql+(RV) :2) ;2 

- rvp^rv 

- c/r*ws*ws*rr 
+ US*RS1*RV 

+ VS+rS2*K\/ 

+ 3*RS*VS2*R\/ 

+ HS*RS2*V\/ 

*• 24^VRU*QL*RS2+ ( RV ) ; 2 
= -f PMR2+WWR1*R0L*(RV) :2) ;2 
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+ PGR*Wt.Rl*RS*USl*RV 
+ PGR*Wt>Rl*RSi4 \ZS2*RV 
+ S+PGR*WGR1+PS/R*US*KV 
+ pGR*WGR2*RS*US?*RU 
+ pGR*WGR2*rs*\/S2*RV 
+ PGR*WGR2=t=RS*WSR2*R*RW 
+ c*PGR*WGR2*PS/R+WS*RW 
+ pGR*WGR3+RS^VS2*RV 

- c*PGR*WGK3+R3/R*WS*RW 
+ PgR2*WWGR+ROOL*R\/ 

- 2*RS*Us1+RV 

- 3 + S*RS/R=»-US^RV 

+ ?/3*EMU*(^-MR + {RV) ;2) :2 

- PML(*RUdMF*RV 

+ 4/3*EMU2+WMR*(PV) ;2 


RSUS* ( R*RW) : 1 

+ RSVS*(R*RW) :2 
+ 2*S*RSUS*Rif, 

- ( vrl)*RQL*(R*RW) ; 2) ;2 

- R\/p*R + RW 

+ US*RR*RWS1 
+ VS+PR+RWS2 
+ US+RS1*R*RW 
+ VS*RS2*H*RW 
+ RS=t:RWSl*RU 
+ RS4.RWS2*KV 
+ RS*HSZ*a*VW 

- VRR+QL+KR2*RWS2 

> VRiJ*QL*RS2=t=(R*RW) ;2 
= pGR*WGKl*RS+uSl=t'R*RW 
+ PGR+WGR1:«RS*VS2*R*RW 

+ g+pgr*wgri*rsus*rw 

- s*pGR*WGR 2*P3 + ^'S*RU 
_ r*P&R*WGR2*RS*WS*RV 
+ s+PGR*WGR 2*FSUS*RW 
+ PGR*WGR3tRS*RWSl=^RU 
+ pgR*WGR 3*RS*PWS2*RV 
+ S*pGR*WGR3*R3US*RW 
~ S*PGR*WGR3*PS*WS+RU 

- r*PGR*WGR3*R3l=^S*RV 
+ PGR2*WWGR=f:RG0L+R*RW 

- ?*RS*US1>*-R*RW 

- 2*RS+WS2*R*PW 

+ pmu*{Wmr*{R’»rw) :2) ;2 

- pmU*RUDMF*R*RW 

+ r|v!U2+WriRt:(R + Ri, ) :2 
+ 1 /2*EMUT*RWS2+(RT) :2 
+ rMUT*RWS22+RT 


MgE is 
POOR QUALErx? 



SIXTH PASS 


14 

RSUS+(HU) :i 

+ rsvs*(Hu) 

- 2 + S*RS/R*l«IS*HW 

~ ( vuH*Ri3L*(iiu) :2) :a 

- RV 2 +HU 

+ US+PU*(HS) :i 
+ VS+RU*(HS) :2 
+ US*RH+US1 
+ VS+RH+US2 

- S/R + fciiS*WS*KH 
+ RS4=UV*(HS) 12 
4- RS+USl+HU 

+ RS*US2*HVJ 

- VRU*0L*RU2*HSJ2 

- - \/RH*QL*RH2*USJ2 

= GMOMS*UU+DPOX 

+ gmoms*uv*qpuy 

+ PGH*WGHl!4.RS*USJL*HlJ 
+ PGH + WGHl*RS*\/S2*HU 
+ S4 pgh*wghi*rs/r*us*hu 

+ PGH*WGH2*RS*US1*HU 
+ PGH*WGH2=*=RS*WSR1 + R*HU 
+ S*PGH*wGH2+RS/R*WS*HvJ 
+ PGH*WGH3*RS*US1*HU 
+ PGH*wGH34^RS*US2*HVJ 

- s*P6H*wGH3*RS/R*wS*hW 
+ PGH2*vJWGH*RQ0L4:HU 

+ Gr0OMS*Ei'1U*US?*UUJ2 
+ 1/2+GM0MS*LMU*RWS2*UWRJ2 
+ 1/2*GM0 HS*SFiU*US 2*VVJ2 
+ l/2+GM0ns*EMU*v4SR2+RUWJ2 
+ GrtOMS+EMUT/CP*US2*USJ2*HU 
+ GF'0MS*EMUT/CP=4WSR2*RtJSJ2*HU 
+ (CAP*(HU) ;2) ;2 

- CAP*HUDFiF*HU 

- (WMT*CAP*(HU) :2) :2 

+ CAYT/CP/CP*(HU*(HS) ;2) :2 

- CAYT/CP/CP*HS2*Wi' 1T*(HU) 52 
+ EMU*< Wf'iH+CHU) :2) :2 

- Er‘lU*HUDMF*HU 

+ ErU2*WMH+(HU) :2 
+ i/ 2+EMUT/CP*US2*{ HH) :2 
+ EPiUT/CP*US22*HH 


RSUS+IR+HW) ;i 

+ RSVS+ ( R*HU ) : 2 

- (\/UH+RQL*{R+HW) ;2) :2 

- R\/2+R*HW 

+ h+us+kw* { hs ) : 1 
+ R*VS*HW*{HS) *.2 
+ US:^RH*RWSl 
+ \/S:t:HH*RWS2 
+ rs:4:r*vw*(hs) :2 


ORIGINAI} PAGP 7<i 

0^- POOE <Sf 
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+ RS+RWSl+HU 
+ RS+RWS2*HVJ 
•» V RU*QL*RRw2*HS J2 
- VRH*QL*RH2*RWSJ2 
= G|VtO«S>t:R*UW*DPOX 
+ G^':OMS*R♦\/W + DPDY 
+ PGH*WGH14.RS*US1*R*HW 
+ PGH+WGHi+Rb*VS2*K*HW 
+ s*pgh*wghi*rsus+hw 

~ S + PGH+U'GH2*RS*WS*HU 

- C*PC>H+WGH2 + RS*WS+HVJ 
+ S+PGH*WGH2*RSUS+HW 

+ PGH*wGH3*kS*RWSi*HU 
+ PGH*WGH3*KS*RWS2*HVJ 
+ S*PGH*WGHi*RSUS*HW 

- S*PGH*WGH3*RS*WS*HU 

- C*P&H*WGH5+'RS*WS*H\/J 
+ PGH2*WWGH*RQ0L*R*HW 
+ gM0MS*EMU*US2+HUWJ2 

+ l/2*GM0MS*Ef1U*RWS2*WWJ2 
+ 1/2+GM0MS*EMU*WSR2+R*R*VVJ2 
+ l/2*GM0MS*EMU=*^tJSR2*RWRJ2 
+ GP>0I“IS*EMUT/CP+US2*USJ2*R>«=HW 
+ GM0MS*EMU7 /CP*WSR2*RWSJ2*R*HW 
+ (C/\P*(R + HW) :2) ;2 

- CAP*HU0MF*R*HW 

- (WMt+CAP*(R*HW) :2) :2 

+ CAYT/CP/CP* ( R*HW+ ( Hs) : 2) ; 2 

- CAYT/CP/CP*HS2*WMT*(R*HW) :? 

+ E[«tu* { WFIH* ( R + HVJ ) : 2 ) : 2 

- E[i''U=fHUDMF*R*HW 

•f LMU2+WMH* < Ri^HV! ) : 2 
+ 1/2*EMUT/CP*RWS2*(HH) J 2 
+ EMUT/CP+RWS22+HH 


BEAM = GYCSZ 


ORIGINMl PXGSTi 
OF POOR QUAlir^ 
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APPENDIX E.- TENSR NOTATION EXAMPLES 


; is used as a separator in these examples. 

Common Notation TENSR Notation 

Scalars 

A ; (|) ; t2 ; V A ; PHI ; TIMEP ; NU 


Covariant Vectors 




A‘I ; 

PHI’R ; V'A 



Contravariant Vectors 


1 1 


A"I ; 

PHI”R ; V"A 



Various Tensors 


A^J 


A"IJ 

or A"I"J 

^ij 


A'lJ 

or A'I'J 

A^ ■ A ^ 


A”I’J 

; A'I"J 

ijm 

Icm 


CAPSIG"IJM’KM 



Covariant Derivative 


♦>i 


PHI, I 




Partial Derivative 


3{|) 

3x^ 


PHI: I 
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Common Notation 


TENSR Notation 


Metric (or Fundamental) Tensor 




@’IJ ; @"IJ 


Kronecker -Delta 


6 ^ 

J 






&’IJ ; &"IJ 


Covariant Derivative with Index Raised 
(or "Contravariant Derivative”) 


im, 

S <i>, 


m 




= PHI! I 
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APPENDIX P,- GLOSSARY 


This glossary includes : 

1) The ;names of the variables appearing in the tensor 
equations of Appendix C and in the final equations 
of Appendix D. 

2) The names of all the subprograms of the GYC s.ystem. 

3) The Fortran names of all the input variables. 

4) The Fortran names of some of the important variables 
of the programs . 

5) The labels used on the program output. 

Many of the entries fit in several of the above categories. 

In the list, several of the names end with . This is 

used to indicate that the name may appear with any of the 

following appearing in place of ; M, Z, P, I, J, K, 1, 2, 

J2, 22, or none of them. The letters refer to the six-point 
finite-difference mesh used in the main solution process as 
indicated in the sketch. 
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The values denoted by names ending in M, Z, and P, at the 
old X , are known. Those ending in I, J, and K may be known 
or unknown depending on the stage of the computation. The 
suffix 1 indicates the first derivative with respect to x 

between Z and J, 2 the first derivative with respect to y 

at Zj J2 the first derivative with respect to y at J, and 

22 the second derivative with respect to y at Z. A name 

that has a different meaning without any such suffix is listed 

both with and without , s-g-j S . On the other hand, a 

variable that carries the same basic meaning without a suffix 
is listed once. For example, HH appears as such in Appendices 
C and D, and as an input variable, as well as in the form HHM , 
HHZ , etc. In any case it signifies h’h’ , so it appears in 
the list only as HH . 

The column labeled Dimensions gives the dimensions of the 
quantity in terms of a density, D , a velocity, V , a length, 
L , a temperature, 0 , and an enthalpy per unit mass, H . 

Q 

(If mechanical units are used for heat, then H = V . The sub- 
stitutions D = ML“3 and V = LT“^ will transform this set 
into a more familiar one.) Absence of any indication means the 
quantity is nondimensional . The symbol W means the elements 
of the vector have the same dimensions as the corresponding 
elements of the vector XV as given in Table 5 . The combination 
W,N used with EPSBV and EPSTV means W if the element is input 
as a positive number and nondimensional if it is input as a 
negative number. Finally, X indicates some other mixture of 
dimensions for different elements. 

At the end of several of the definitions there are expres- 
sions in square brackets. These are evaluations in terms of 
Mach number, Reynolds number, etc. of the variables when they 
are interpreted as nondimensional quantities under the following 
conditions : 
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1) The nondimensionalization is done with respect to a 

length, , and with respect to a reference density, 

pp , velocity, , temperature, T^ , and enthalpy 

per unit mass, h^ , that all refer to the same con- 
dition or point in the flow. 

2) The reference state for the transport coefficients is 
that same condition, i.e., T^f = T^ . 

An example of how they are derived is given for the case of 
GMOMS in Section 2 . 

Notes to the Glossary appear at the end. 


Name Dimensions 


Definition 


A 

AHH 

AHU 

AMAT 

ARU 

AV_ 

B 

BDGP 

BDV 

BETA 

BHH 

BHU 


Modeling parameter, a . See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. 

X Matrix A (Section 3). 

Modeling parameter. See Note 2. 

W Vectors of the unknowns at the new x . 

AVM, AVZ, and AVP would more logically 
be called XVI, XVJ, and XVK. 

Modeling parameter, b . See Note 2. 

1/(Dv2) -BDV/(yp) 

Modeling parameter. See Note 2. 

Modeling parameter, B . See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. ,■ 
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Name 

BLAM_ 

BLAMO 

BMAT 

BREAK 

POINTS 

BRU 

• BUFAC 
C 

CAP_ 

CAY_ 

CAYRF 

CAYT_ 

CP 

CHANGE 

STEP 

CLAMB 

CMAT 

CMNTS 

CONGL 

CP 


Dimensions Definition 

L Turbulence scale length. See Tables 3 

and 5 and Note 1. 

L Scale parameter. See Note 1. 

X Matrix B (Section 3)- 

L TBRKV 

Modeling parameter. See Note 2. 

Maximum ratio by which EPSTV criteria 
may be exceeded without causing- backup 

DIMPLS cos 0^, 

DVL k/Cp 

DHVL/0 Heat conductivity, k . 

DHVL/0 Reference value of k : its value at 

TMRP. See Note 3. 

[1/ (RerPr ) ] 

p 

DHVL/0^ Derivative of k with respect to 
temperature . 

Skin friction coefficient. 


PER W,N EPSTV 


Scale parameter. See Note 1. 

X Matrix C (Section 3)* 

Up to 76 characters of comments to be 
printed on the' title pages. 

Cone half-angle, 0^ , in degrees. 
See Note 5- 

H/0 Specific heat at constant pressure, 

Cp , assumed constant. 

[1.3 
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Name 


Dimensions 


Definition 


CSUTH 

9 

Constant in the viscosity law. See 
Mote 3. 

DELT 

L 

Step size in the x direction. Ax . 

DELTA* 

L 

Displacement thickness. 

DELTI 

L 

Input value of DELT. See Note 7. 

DIMPL 


Important; see Note 5 . 

Dimension flag. A value of zero signi- 
fies flat-plate mode; a value of one 
signifies axisymmetric mode. Any other 
input value is converted to one. 

DIN 


Scale parameter. See Note 1. 

DH/DY 

H/L 

The wall value of ^ . 

9y 

DLT99 

L 

Value of y for which u is 99 percent 
of its free stream value. 

DOUT - 


Scale parameter. See Note 1. 

DPDX 

Dv2/L 

Partial derivative of p with respect 
to X . 

DPDY 

DV^/L 

Partial derivative of p with respect 
to y . 

DSM 

L 

SZ - SM = h (Section 3). 

DSP 

L 

SP - SZ = h^ (Section 3). 

DST 

L 

SP - SM = h.[. (Section 3)* 

DTPMX 


Maximum ratio by which Ax may increase 
from one step to the next. 

DTI* 

L 

Velocity thickness. 

DTMAX 

L 

Maximum Ax . 

DTMIN 

L 

Minimum Ax . 



Name 


Dimensions 


Definition 


DU/DX V/L 

DVEC X 

DZPMN 

DZPMX 

DZRMX 

ECMN 


ECMNI 

ECMX 

ECMXI 


EDGE 

TOLERANCE 

W,N 

EDGEV 

W 

EMS 

DVL 

EMST 

DVL/e 

EMU_ 

DVL 

EMURF 

DVL 

EMUT_ 

DVL/0 

EPSBV 

W,N 


Vector D (Section 3)* 

DZEMN^TYPL is the minimum Ay . 

DZPMX*TYPL is the maximum Ay . 

Maximum ratio of two adjacent, values 
of Ay . 

Governs dropping of points from the 
profile . 

Initial value of ECMN at each step. 
Governs adding points in the profile. 
Initial value of ECMXI at each step . 
EPSBV 


Free stream values of dependent 
variables . 

Second coefficient of viscosity, y* . 

Derivative of y* with respect to 
temperature . 

First coefficient of viscosity, y . 

Reference value of y : its value at 

TMRP. See Note 3. 

[l/Re^f.] 

Derivative of y with respect to 
temperature . 

The solution is extended in the radial 
direction until the difference between 
the calculated value of each parabolic 
variable and its asymptotic value .is 
less than EPSBV(I) or -EPSBV(I)^fGMAXV(I) , 
whichever is positive . 
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Name 


Dimensions 


Definition 


EPSNV 


EPaw 


EV 


W When ABS(GMAXV(I) ) is less than 

ABS(EPSNV(I) ) , the corresponding 
variable is not checked in the deter- 
mination of DELTj nor in considering 
where to insert or drop points, nor 
in checking the outer boundary con- 
ditions if EPSBV(I) is negative. 
Furthermore, in the determination of 
Ax , if EPSNV(I) is negative, the 
change in the corresponding variable 
is compared with the larger of 
-EPSNV(I) and the criterion derived 
from EPSTV(I) . 

¥,N Ax is controlled by attempting to 

keep the change, TCMXVCI), in each of 
the parabolic variables (in the axial 
direction) below EPSTV(I) or -EPSTV* 
GMAXV(I), whichever is positive. 

See Table 2, 


EXPMU 


Exponent in the viscosity law. See 
Note 3. 


FAC 9 9 


See Note 1. 


FCMX 

FCURI 

FHM 

FHP 

FHZ 

FINTV 


Maximum factor by which ECMNI and 
ECMXI may be increased through repeated 
applications of FCURI. 

Factor by which ECMN and ECMX are 
Increased when JMAX impedes adding 
points in the profile. 

1/L -H_ (Section 3)* 

1/L (Section 3)* 

1/L H_ - (Section 3). 

L Spreads of the dependent variables 

based on their Integrals (Output label: 

INT SPREAD) . 
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Name 


Dimensions 


Definition 


FMAXV 

FUP 

FUPV 


¥ Maxima over y of the dependent 

variables (output label: MAXIMUM). 

Upwind differencing parameter. See 
Note 4. 

Upwind differencing parameter. See 
Note 4. 


FV 


See Table 3. 


FZMXV 


L Locations of the maxima (output label: 

Y AT MAX) . 


GLOBAL MAX W GMAXV 


GMAT 


X Matrix r (Section 3)* 


GMAXV 

GMOMS 


GV 


W Maxima over y and x of the dependent 

variables (output label: GLOBAL MAX). 

Important: see Note 6. 

H/V^ The reciprocal of the mechanical 
equivalent of heat . 

(y-DmJ 

See Table 4. 


GYCAA 

GYCAL 

GYCAQ 

GYCBI 

GYCBO 


Subroutine which adjusts the spacing of 
points in the y direction. 

Subroutine which controls output and 
auxiliary computations. 

Subroutine which computes spreads and 
the next step size and which initializes 
for the next step. 

Subroutine which implements the wall 
boundary conditions. 

Subroutine which implements the outer 
boundary conditions. 
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Name 

GYCCM 

GYCCP 

GYCDZ 

GYCEW 

GYCFP 

6YCIC 

GYCIN 

GYCIP 

GYCLI 

GYCMC 

GYCMI 

6YCMP 

GYCNS 


Dimensions 


Definition 


COMMON file. 

Subroutine which calculates the pressure 
distribution and the density distri- 
bution when they need to be determined 
separately. 

Subroutine which determines the position 
in y of a new point which is to be 
inserted . 

Subroutine which interpolates to deter- 
mine values of the specified free-stream 
variables at the current x . 

Subroutine which transfers from STOR, 
or calculates, the values in the section 
of COMMON from QP to USKl. 

Subroutine which handles some of the 
initialization for a new run. 

Subroutine which controls start or 
restart procedures. 

Subroutine which determines values of 
the dependent variables at a point 
inserted in the profile. 

Subroutine for linear interpolation. 

Subroutine which controls the compu- 
tation of AMAT, BMAT, CMAT, and DVEC . 

Subroutine which inverts 2 • 2 and 
3 • 3 matrices. 

Subroutine which computes a matrix 
product and optionally subtracts the 
product from a third matrix. 

Subroutine which flags certain unlikely 
errors. 
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Name 

GYCOI 

GYCON 

GYCOT 

GYCPP 

' GYCPI 

GYCPS 

GYCPI 

GYCP2 

GYCP3 

GYCP4 

GYCP5 

GYCP6 

GYCQD 

GYCRI 

GYCRP 

GYCRV 


Dimensions Definition 

Subroutine which produces a minor 
printout and the first part of a 
major printout. 

Main program of the GYC system. 

j Subroutine which produces the three 

profile tables of a major printout 
and also prints any error messages. 

Plcticious subroutine which would 
write to a file for later plotting 
of profiles. 

Subroutine which prints the first two 
pages of a new or restarted run. 

Ficticious subroutine which would set 
up plot files. 


Subroutine 

to 

GYCMC 

for 

the 

second pass. 

Subroutine 

to 

GYCMC 

for 

the 

first 

pass . 

Subroutine 

to 

GYCMC 

for 

the 

third 

pass . 

Subroutine 

to 

GYCMC 

for 

the 

fourth pass . 

Subroutine 

to 

GYCMC 

for 

the 

fifth 

pass . 

Subroutine 

to 

GYCMC 

for 

the 

sixth 

pass. 


Subroutine which calculates various 
integrals and determines maxima. 


Subroutine which reads all but the first 
input card. 

Ficticious subroutine which would write 
to a file for later plotting of "running 
variables" as a function of x . 

Subroutine which determines the normal 
component of the mean velocity. 
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Name 


Dimensions 


Definition 


GYCSC 

L 

Function which determines the value of 
the turbulent scale. See Note 1. 

GYCSS 


Subroutine which controls the main 
solution . process . 

GYCSZ 

L 

GYCSC(SZ) 

GYCTB 


Subroutine which makes a final adjust- 
ment to the step size. Ax , if needed 
to hit a break point. 

GYCTC 


Subroutine which evaluates k and y 
and their derivatives with respect to 
temperature. See Note 3- 

H 

H 

H . See Tables 3 and 5. 

HE 

H 

Free-stream value of h . See Table 2 

HG 

H/L 

hy at the wall. 

HH_ 


h’h' . See Tables 4 and 5* 

HHDMP 

1/l2 

(AHH + BHH»REL)/BLAMZ/BLAMZ 

HS_ 

H 

h . See Table 5* 

HSTAG 

H 

The stagnation enthalpy is ABS(HSTAG). 
A positive value indicates that the 
wall value of the mean enthalpy is 
specified on Cards l6.n. A negative 
value indicates that the wall value of 
the gradient of the mean enthalpy is 
specified on Cards l6.n. 

[i + 1=1 mI] 

HSW 

H,H/L 

The wall value of ii or of . See 

Table 2. ^ 

HU_ 

HV 

h'u’ . See Tables 4 and 5* 

HUDMP 

1/l2 

(AHU + BHU*REL)/BLAMZ/BLAMZ 

HV_ 

HV 

1 

h'v' . See Tables 4 and 5- 


80 



Name 


Dimensions 


Definition 


HW_ 

I 

INPLG 

INT SPREAD 

lOFAC 

ITEST 

J 

JMAX 

JMXMX 

JTOP 

LAST DX 
LBURF 

LDELT 

LDETZ 

LFOLP 


HV h'w' . See Tables 4 and 5- 

Index almost always used to range over 
the dependent variables. 

1 flags a new run (all cards are read); 

2 flags a restart with changes (Cards 1 
through 29 are read); 0 flags end of 
job. See Note 7* 

L ' FINTV 

Scale parameter. See Note 1. 

Index of variable to which PCPMX and 
PCPSP apply. 

Index used to range across the profile. 
J=1 corresponds to the wall; J=JTOP 
corresponds to the edge of the boundary 
layer (free stream) . 

Maximum number of points to be allowed 
across the profile, i.e., the upper 
limit on JTOP. 

Maximum permissible value of JMAX. 

Current number of points across the 
profile . 

L Ax for the step just completed. 

True indicates that the step is to be 
tried again with a smaller Ax . 

True Indicates that Ax is less than 
its minimum. 

True indicates that GYCMI has tried to 
invert a singular, or nearly singular, 
matrix. 

True indicates that full output is to 
occur. 


8l 



Name 


Dimensions 


Definition 


LPOPP 

LHMWK 

LIOLP 

LIOPP 

LJMXA 

LNWPG 

LPKRF 

LPOSV 


True indicates that output to a profile 
plot file would occur. 

True indicates that the wall value of 
enthalpy is specified; false that its 
gradient is specified. 

True indicates that a minor printout is 
to occur. 

True indicates that output to a running 
plot file would occur. 

True indicates that more than JMAX 
points are needed to satisfy the curva- 
ture requirements. 

True indicates that the next, printout 
should start on a new page . 

True indicates that more than JMAX 
points are needed to satisfy the outer 
boundary conditions. 

True Indicates a non-negative dependent 
variable. 


LPPPF 


True would indicate that a profile file 
was full. 


LRHOF True indicates -a problem in the calcu- 
lation of p : either h*^ - h^ or 

too many iterations in GYCCP. 


LRORG 

LRPPP 

LRSTP 


'True Indicates that points have been 
added or deleted in the current step. 

True would indicate that a running file 
was full. 

True indicates that a restart is in 
progress . 
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Name 


Dimensions 


Definition 


LSTFL 


True indicates that a new run is 
starting. 

LTRNP 


True indicates that the^ run is termi- 
nating . 

LZHPF 


True indicates that y has reached 
YMAX. 

LZHRP 


True indicates that correlations 
involving p’ and h’ are to he set 



to zero. 

MAX CHANGE 

W 

TCMXV 

MAXIMUM 

W 

FMAXV 

NEXT DX 

L 

Ax for the next step. 

NPOLP 


Steps between major printouts. 

NPOPP 


Steps between fictitious writes to a 
profile plot file. 

NIOLP 


Steps between minor printouts. 

NIOPP 


Steps between fictitious writes to a 
running plot file. 

NOISE 

W 

EPSNV 

THRESHOLD 



NPASS 


Number of the pass through the procedure 
of GYCSS. 

NRUN 


Run number. 

NSKIP 


BLYAA is skipped for the first NSKIP 
steps whether starting or restarting. 

NSTEP 


Number of the current step in the x 
direction. 

NSTMX 


The run is halted if the number of steps 
reaches NSTMX. See Note J. 
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, Name 


Dimensions 


Definition 


NSVRT 

NVARP 

NVART 

NWQRD 

NWSAV 

NWSGM 

NWSRD 

NWVEC 

OPAC 

OMBET 

OUTPUT 

CONTROLS 

P 

PCPMX 

PCPSP 

PCTMV 


NVART’ plus one for the dependent vari- 
able y (SJ; 22. 

Number of parabolic variables; 17. 

Total number of dependent variables 
saved for the old and the new value of 

X ; 21. ' ■ 

Number of words in COMMON representing 
the M-Ij the Z-J, or the P-K conditions; 
59. 

Total number of variables saved; 43. 
NWSAV plus the size of GMAT; 52. 

NWSAV; 43. 

Number of unknowns in a pass; 2 in the 
first pass, 3 in the others. 

Scale parameter. See Note 1. 

1 - BETA. 

Parameters input on Card 17. 


DV^ p . See Note 6. 

If PMAXV(ITEST) has changed more than 
PCPMX times the value it had the last 
time both BLYOT and BLYPP were called, 
a call to each is forced. 

If ZPSM(ITEST) has changed more than 
PCPSP times the value it had the last 
time both BLYOT and BLYPP were called, 
a call to each is forced. 

ZPSM(I) is the largest value of y 
for which XV(I) = PCTMV(I)5^PMAXV(I) . 
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Name 


Dimensions 


Def-inition 


PCT SPREAD 


L 

ZPSM 



PE 


DV^ 

Free-stream value of 
and Note J. 

p . See Table 2 

PGH 



Modeling parameter. 

See Note 

2. 

PGH2 , 



Modeling parameter. 

See Note 

2 . 

'PGR 



Modeling parameter. 

See Note 

2. 

PGR2 



Modeling parameter. 

See Note 

,2. 

PGU 



Modeling parameter. 

See Note 

2. 

PGU2 



Modeling parameter. 

See Note 

2. 

PMH 



Modeling parameter. 

See Note 

2. 

PMH2 



Modeling parameter. 

See Note 

2. 

PMR 



Modeling parameter. 

See Note 

2. 

PMR2 



Modeling parameter. 

See Note 

2. 

PMU 



Modeling parameter. 

See Note 

2. 

PMU2 



Modeling parameter. 

See Note 

2. 

PS_ 

Dv2 


p . See Note 6. 



PSTAG 

DV2 


Important: see Note 6. 

The stagnation pressure is ABS(PSTAG). 
A positive value Indicates a normal 


run; a negative value indicates that 
all correlations involving p’ and h’ 
are zeroed out at each step. (This is 
intended as a debugging device.) 

(i + 

PTHM Modeling parameter. See Note 2. 

PTS JTOP 
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Name 

Dimensions 

Definition 

PTUM 


Modeling parameter. See Note 2. 

i ' 

Q_ 

V 

q = ~wu'u' + v’v' + w’w' 

QL_ 

W 

qA 

QQ 

V2 

UUZ + WZ + WWZ 

QQJ . 

V2 

UUJ + WJ + WWJ 

QTOT 

DHV 

Total heat transfer rate in the y 
direction. 

QWALL 

DHV 

Heat transfer rate out of the wall. 

R_ 

L 

Distance to the axis of the cone. 
See Note 5- 

RBLS_ 

Dl2 

pa2 

RDT 

1/L 

1/Ax 

REL 


pqA/y 

RETHETA 


Reynolds number based on momentum 
thickness . 

REXW 


Reynolds number based on wetted length. 

RGAS 

v2/0 

Important: see Note 6. 

Gas constant . 

i/(ym2) 

RH_ 

DH 

p’h’ 

RHO 

D 

P 

RNOSE 

L 

Origin radius, r^ . See Note 5* 

ROTAT 

V/L 

Important: see Note 5- 

Rate of rotation of the cone in radians 

per unit time. 
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Name 

Dimensions 

Definition 

RQL_ 

DVL 

pqA 

RQOL 

DV/L 

pq/A 

RR_ 

d2 

P’P’ 

RRW2 

D^V/L 

(PP’W )y 

RS_ 

D 

P 

RSUS 

DV 

RS*US 

RSVS 

DV 

RS«VS 

RT_ 

D0 

ptipt 

RTW2 

DV0/L 

(pT’w"* )y 

RU__ 

DV 

p'u' . See Tables and 5* 

RUDMP 

1/l2 

(ARU + BRUsRED/BLAMZ/BLAMZ 

RUN 


NRUN 

RUWJ2 


(Ru'w')y at the new value of x . 

RV_ 

DV 

p’v’ . See Tables 4 and 5. 

RWJ2 

v2 

(Rv'w')y at the new value of x . 

RW_ 

DV 

p'w’ . See Tables 4 and 5- 

RWRJ2 

V^L 

(R'^w'wMy at the new value of x . 

RWS_ 

DV 

Rw 

S 

■ 

DIMPLE sin 0^ 


L 

Normal independent variable ^ y . 

SFVPL 


Subroutine which fills a vector. 

SPVMV 


Subroutine which moves the contents 
one vector to another. 
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Name 


Dimensions 


Definition 


SHAPE 


Displacement thickness divided by 
momentum thickness. ‘ 

SHM 

1/l2 

. 1/ ^h_h^^ (Section 3). 

SHP 

1/L^ 

l/(h+h^^ (Section 3)- 

SPREAD 

PARAMETER 


PCTMV 

SPRO 

L 

ZPSM(IOPAC). See Note 1. 

STANT 


Stanton number. 

STEP 


NSTEP 

STOR 

X 

Storage array mimicking an external 
device . 

TAULAM 

DV^ 

Laminar stress. 

TAUTOT 

DV^ 

Total stress. 

TAUW 

Dv2 

Wall stress. 

TBRKV 

L 

Time break vector. When TIME reaches 
the absolute value of one of the ele- 
ments of TBRKV, calls to GYCOT and 
GYCPF are forced. If the element is 
negative, the run is halted. See 
Note 7. 

TCMXV 

W 

Maxima over y of the changes in the 
dependent variables during the last 
step (output label: MAX CHANGE). 

THETA 

L 

Momentum thickness. 

TIME 

L 

Streamwise' independent variable, x . 
See Note 7 • 

TMRF 

0 

Reference temperature for y and k 
See Note 3. 

[1.] 
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Name 


Dimensions 


Definition 


TT 

TU_ 

TV_ 

TYPL 

U 

UE 

US_ 

uu_ 

uv_ 

uw_ 

UWRJ2 

V 

VARNM 

VBS 

VHH 

VRH 

VRR 

VRU 

VS_ 

VSW 

VUH 


,02 

V0 

V0 

V0 

L Typical length. See Note 1. 

V u . See Tables 3 and 5- 

V Free-stream value of u . See Table 2. 

V u . See Table 5* 

u'u' . See Tables 3 and 5* 

y2 u"* v’’’ . See Tables 3 and 5. 

u’w’ . See Tables 4 and 5- 

V^/L^ (u’w'/R)„ at the new value of x . 

«y 

V V . See Table 3. 

Vector of variable names. 

Modeling parameter. See Note 2. 
Modeling parameter. See Note 2. 
Modeling parameter. See Note 2. 
Modeling parameter. See Note 2. 
Modeling parameter. See Note 2. 

V V 

V The wall value of v . See Table 2 . 
Modeling parameter. See Note 2. 
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Name 


Dimensions 


Definition 


VUU 

vw_ 

VWRJ2 

W 

WGHl 

¥GH2 

WGH3 

WGRl 

WGR2 

WGR3 

¥GU 

¥MD 

¥MH 

¥MR 

¥MT 

¥S_ 

¥SR_ 

¥¥_ 

¥¥D1 

W¥D2 


Modeling parameter. See Note 2. 
v' v' . See Tables 3 and 5- 

v2 v'w* . See Tables 4 and 5* 

V^/L^ (v’w'/R)„ at the new value of x . 

V 

V w . See Tables 3 and 5- 

Modeling parameter computed from 
constraint. See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. 

Modeling' parameter computed from 
constraint. See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter* See Note 2. 

Modeling parameter. See Note 2. 

Modeling parameter. See Note 2. 

V w . See Table 5. 

V/L w/R 

w'w' . See Tables 3 and 5* 

Modeling parameter. See Note 2. 
Modeling parameter. See Note 2. 
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Name 

WWGH 

WWGR 

WWGU 

WWHl 

WWH2 

WWRl 

WWR2 

WWUl 

WWU2 

X 

XSTRT 

XV_ 

xw 

XWET 

Y 

Y AT MX 
YMX 

ZPSM 


Dimensions Definition 


Modeling 

parameter. 

See 

Note 

2 . 

Modeling 

parameter . 

See 

Note 

2. 

Modeling 

parameter . 

See 

Note 

2. 

Modeling 

parameter. 

See 

Note 

2. 

Modeling 

parameter . 

See 

Note 

2. 

Modeling 

parameter. 

See 

Note 

2. 

Modeling 

parameter. 

See 

Note 

2. 

Modeling 

parameter. 

See 

Note 

2. 

Modeling 

parameter. 

See 

Note 

2. 

X . See 

Table 2 . 





L Initial value of XWET. See Note 7- 

W Vectors of the unknowns at the old x 

Compare AV, See Table 5* 

L Wetted length. 

L Wetted length (output label: XW) . 

L y . See Tables 3 and 4. 

L PZMXV 

L The maximum value of y to which the 

solution may spread. 

L Spreads of the dependent variables 

based on PCTMV' (output label: PCT 

SPREAD) . 
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Note 1. 


The turbulence scale A that enters in many of the models 
is computed aS' follows. First, a typical length is computed 
from 

TYPL = PAC99«DLT99 + OFAC*ZPSM(IOFAC) 

where DLT99 is the value of y where u is 99 percent of the 
free-stream value; ZPSM is the vector of spreads of the depend- 
ent variables based on PCTMV; and PAC99j OFAC, and lOPAC are 
inputs. (ZPSM(IOPAC) is included' in the output with the label 
SPRO. ) 

The value of A is made equal to BLAMO at the wall. This 
parameter is ordinarily zero but may be made positive to simu- 
late roughness. Off the wall, A increases with slope DIN 
until it reaches a maximum of CLAMB»TYPL. For y > DLT99> A 
increases again with slope DOUT. This parameter is ordinarily 
■zero but should be positive if there is free-stream turbulence. 

All the parameters determining A are nondimensional 
except BLAMO. 

Note 2. 

* 

The modeling parameters are defined by their appearance in 
the models, as given in Appendix B. They are all nondimensional. 
Two of them, WGHl and WGRl, are not input; instead, they are 
computed from the constraints given at the end of Appendix B. 
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Note 3- 


The viscosity is determined from 


V 




rf T 


1 + C/T 


rf . 


where iij,f = EMURP , T^f = TMRF , w = EXPMU , and C = CSUTH . 
With 0 ) = 1/2 , this is the Sutherland relation. With C = 0 ^ 
it is the power-law relation. With both w = 0 and C = 0 , 
it gives constant jl . 

The heat conductivity is determined from 


k 



y 

Urf 


Since c is constant, this is equivalent to assuming that the 

’ Jr ' 

Prandtl number is constant. 


Note 

Upwind differencing is suppressed when PUP = 0 , used 
when PUP = 1 . In this case, the coefficient, pv , is 
evaluated as (1 - PUPV) times the central-point value and FUPV 
times the upwind-point value. If PUP = 0 , the value of 
FUPV is immaterial. Values of PUP other than 0 and 1 produce 
a mixture of standard and upwind differencing. 


Note 5 « 

. If the dimension flag, DIMPL, is input as zero, the program 
operates in flat-plate mode. In this case, RNOSE and C0N6L are 
not used and the following acquire different meanings : 

R_ (Dimensionless) is the number 1. 

ROTAT (Dimension: V) is the sideslip velocity of the surface. 
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Thus, the value of w at the surface is 


ROTATSR 

(where R is evaluated at the surface) for either value of 
DIMFL. 

Note 6 . 

If GMOMS is input as zero, the program operates in Mach 
number zero mode. In this case, the following acquire different 
meanings : 

GMOMS (Dimensionless) is the number 0. The units for 
expressing heat quantities may still be chosen at will. 

A conversion factor between mechanical and heat units is 
not needed in this case. 

P, PE, and PS_ (Dimensions: DV^) all represent pressure 
differences instead of absolute pressures. The datum for 
zero values for the outputs is whatever is used for the 
input (PE). If the free-stream value is constant, it 
would ordinarily be input as zero, and the output pressures 
would represent departures therefrom. If PE is input as 
the absolute free-stream pressure, the output pressures 
are absolute too. 

PSTAG (Dimensionless) is one of the numbers ±1 . This is 
used only to flag runs for which correlations involving 
p'’ and h’ are to be zeroed out, 

RGAS (Dimension: D) is the density corresponding to HSTAG. 


9 ^ 



Note 7 » 

On restart (INPLG = 2) two inputs, DELTI and XSTRT, are 
ignored if input as zero; the current values of DELT and XWET 
in the halted run are used Instead, On the other hand, the 
input values will be used if they are not zero. The input 
value of TIME (x) is always used, so if it is desired to have 
it run continuously the final value from the halted run must 
be supplied as input. The time break values, TBRKV, are 
interpreted in terms of TIME, so they may need to be altered 
on restart. If the run was halted by the use of a negative 
time break, that should probably be changed for restart. If 
the run was halted by NSTMX, that should be changed for restart. 
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Table 1. Input Cards 


Card 

1 

INFLG 

, NSKIP 


(I1,3X,I4) 

Card 

2 

CMNTS 

, NRUN 


(19A4,i4) 

Card 

3 

DIMPL 

, RNOSE, CONGL, ROTAT, TIME 


(5P8,0) 

Card 

4 

EGAS, 

CP, HSTAG, PSTAG, GMOMS, XSTRT 


(6P8.0) 

Card 

5 

TMRP, 

EMURP, CAYRP, CSUTH, EXPMU 


(5P8. O) 

Card 

6 

CLAMB 

, DIN, DOUT, BLAMO, PAC99, OPAC, 

I OPAC 






(6p8.0,i4) 

Card 

7 

A, B, 

BETA 


(35'8.0) 

Card 

8 

AHH, . 

AHU, ARU, BHH, BHU, BRU 


(6P8.0) 

Card 

9 

vuu, ■ 

VUH, VHH, VRU, VRH, VRR 


(6P8.0) 

Card 

10 

PMU, : 

PMU2, PMH2, PMR2, PTHM, PTUM 


C6P8.0) 

Card 

11 

PGH, 

PGH2, PGR, PGR2, PGU, PGU2 


(6F8.0) 

Card 

12 

, WGU, ■ 

WMH, WMR, WMT 


(4P8.0) 

Card 

13 

WGH2, 

WGH3, WGR2, WGR3 


C4P8.0) 

Card 

14 

WWUl, 

WWRl, WWHl, WWGU, WWGR, WWGH 


C6P8.0) 

Card 

15 

EV 



(6F8.0) 

Cards l6.n 

X, UE 

, HE, PE, VSW, HS¥ 


(6p8.0) 

Card 

17 

NIOLP 

, NIOPP, NFOLP, NPOPP, PCPMX, PCPSP 

, ITEST 





(4l4,2P8.0,l4) 

Card 

18 

TBRKV 


(10F8.0) 

Card 

19 

NSTMX 

, JMAX, YMAX, PUP, PUPV 


(2I4,3P8.0) 

Card 

20 

ECMNI 

, ECMXI, FCURI, PCMX, DZPMN, DZPMX, 

DZRMX 






(7P8.0) 

Card 

21 

DELTI 

, DTMIN, DTMAX, DTFMX, BUPAC 


(5P8.0) 
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Cards 22/23 
Cards 24/25 
Cards 26/27 
Cards 28/29 
Card 30 
Cards 31. j 
Card 32 
Cards 33. j 


Table 1. Input Cards (Continued) 


PCTMV 

EPSTV 

EPSBV 

EPSNV 

EV 

y, U, V, W, H, 
GV 

Y, UW, W, HU, 


UU, W, WW, UV, 
HV, HW, RU, RV, 


(8P8.0/9E8.0) 
(8P8.0/9P8.0) 
(8F8.0/9P8.0) 
(8P8.0/9P8.0) 
(10P8.0) 
BLAM (10P8.0) 

(10P8.0) 
RW, HH (10P8.0) 
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Table 2. Edge Values and Wall Values 


Name 

Multiplied By 

Gives 

Dimensions 

X 

EV(1) 

X 

(L) 

UE 

EV(2) 

^e 

(V) 

HE 

EV(3) 

^e 

(H) 

PE 

EV(4) 

^e 

(DV^) 

VSW 

EV(5) 

% 

(V) 



f ^w 

(H) 

HS¥ 

EV(6) < 


(H/L) 


1 2 

If hg is input as zero^ it is calculated from hg = hg - — u^ 
where hg is the stagnation enthalpy. 

If GMOMS > 0 and Pg is input as zero, it is calculated from 
Pe “ Ps where Pg is the stagnation pressure. 
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Table 3* Initial Conditions ^ First Set 


Name 

Multiplied 

By 

Gives 

Dimensions 

I 

FV(1) 


y 

(L) 

U 

FV(2) 


u 

(V) 

V 

PV(3) 


V 

(V) 

w 

PV(i|) 


w 

(V) 

H 

FV(5) 


h 

(H) 

uu 

FV(6) 


u’u' 

(y2) 

w 

PV(7) 


v'v' 

(V^) 

W¥ 

FV(8) 


w’w' 

(V^) 

UV 

PV(9) 


u'v* 

(V^) 

BLAM 

PV(IO) 


A 

(L) 


Table 4. Initial 

Conditions, Second 

Set 

Name 

Multiplied 

By 

Gives 

Dimensions 

I 

GV(1) 


y 

(L) 

UW . 

GV(2) 


u’w’ 

(V^) 

W 

GV(3) 


v' w' 

(V^) 

HU 

GV(4) 


h'u* 

(HV) 

HV 

GV(5) 


h' V* 

(HV) 

HW 

GV(6) 


h'w' 

(HV) 

RU 

GV(7) 


p 'u' 

(DV) 

RV 

GV(8) 


p ' v' 

(DV) 

RW 

GV(9) 


p 'w* 

(DV) 

HH 

GV(IO) 


h'h' 

(H^) 
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Table 5« Dependent Variables 


Name 

Dimensions 

Output Label 

Meaning 

WS = XV(1) 

(V) 

W 

' w 

VW = XV(2) 

(V^) 

VW 

v’w' 

US = XV(3) 

(V) 

u 

u 

UV = XV(4) 

(v2) 

UV 

u'v’ 

UW = XV(5) 

(V^) 

UW 

u'w’ 

UU. = XV(6) 

(V^) 

UU 


W = XV(7) 

(V2) 

w 

v'v' 

WW = XV(8) 

(V^) 

WW 

w’w' 

HS = XV(9) 

(H) 

H 

h 

HV = XV(IO) 

(HV) 

. HV 

h’v' 

HH = XV(ll) 

(H^) 

HH 

h’h’ 

RU = XV(12) 

(DV) 

RU 

p ’u ' 

RV = XY(13) 

(DV) 

RV 

p 'v' 

RW = XV(14) 

(DV) 

RW 

p ' w' 

HU = XV(15) 

(HV) 

HU 

h'u‘ 

H¥ = XV(16) 

(HV) 

HW 

h’ w’ 

BLAM = XV (17) 

(L) 

BLAM 

A 
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